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Magnetoelastic  Surface  Waves  In  Half-Space  With  Surface  Electric  Current 

917F0082A  IZVESTIYA  AKADEMII  NAUK  ARMYANSKOY  SSR:  SERIYA  MEKHANIKA  in 
Russian  Vol  43  No  1,  1990  (manuscript  received  25  Oct  88)  pp  3-9 

[Article  by  M.  V.  Belubekyan,  Institute  of  Mechanics,  ArmSSR  Academy  of 
Sciences] 

[Text]  The  earliest  research  on  magnetoelastic  surface  waves  in  the 
presence  of  a  constant  magnetic  field  was  done  by  S.  Kaliskiy  [1].  This 
work  has  subsequently  been  continued  by  many  researchers. 

In  contrast  to  other  papers,  this  one  considers  the  case  where  a  magnet¬ 
ic  field  is  set  up  by  electric  current  flowing  over  the  surface  of  a 
half-space.  Studies  of  the  propagation  of  elastic  waves  in  the  presence 
of  surface  current  may  have  application  in  nondestructive  testing  of 
materials  [ 2 ] . 

1.  Let  an  elastic  half -space  occupy  region  ^  0  in  rectangular  coor¬ 
dinate  system  (Xi,  Xa,  Xg).  Electric  current  with  surface  density  Jq 
flows  along  boundary  Xa  =  0  in  the  direction  of  axis  Ox^.  In  the  ini¬ 
tial  (unperturbed)  state,  the  magnetic  fields  caused  by  Jo  will  be 
constant  in  regions  Xa  >  0  and  Xa  <  0,  undergoing  a  discontinuity  (jump) 
on  surface  Xa  =  0.  The  magnetic  field  intensities  must  be  defined  to 
satisfy  boundary  condition 


_  _  A  4"  , 

Here  Ho*’,  Ho  are  magnetic  field  intensity  vectors  where  Xa  <  0,  Xa  >  0 
respectively,  1,  unit  vectors  in  directions  x^,  Xa  and  Xg 

respectively. 

The  following  three  versions  are  possible  for  definition  of  a  magnetic 
field  satisfying  these  conditions: 
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-  A  _  A  (1,1) 

/7o=0.  0  ^  —Jo  O:--’) 

R^-=z2hX  ^^''>=0 


A  model  for  versions  (1.2)  and  (1.3)  is  a  medium  in  the  form  of  a  layer 
of  rather  considerable  thickness  with  current  flowing  over  its  surfaces 
with  surface  density  Jo.  Where  currents  over  the  surfaces  have  the  same 
direction,  the  magnetic  field  cannot  penetrate  into  the  layer,  and  we 
arrive  at  version  (1.2).  Version  (1.3)  is  the  result  in  the  case  where 
surface  currents  are  oppositely  directed. 

According  to  (1.1)-(1.3,  the  components  of  the  Maxwell  tensors  for  the 
corresponding  versions  take  the  form  (only  non-zero  components  are 
given): 


7  _T(<f)_ _ \.H'i  T  — _ !- T  T  <  '= -L  A/- 


(1.4) 

(!.6) 

(1.6) 


Here  and  below,  the  permeabilities  of  the  material  of  the  half-space 
Xs  S  9  and  of  the  medium  outside  of  the  half-space  are  assumed  to  be 
equal  to  unity. 

In  the  first  version  (1.4),  the  components  of  the  Maxwell  tensor  are 
continuous  on  boundary  Xa  =  0,  and  therefore  the  initial  stressed  state 
is  zero  =0).  In  the  second  (1.5)  and  third  (1.6)  versions,  compo¬ 
nent  T330  undergoes  a  discontinuity  on  boundary  Xa  =  0,  giving  rise  to 
initial  stresses. 


To  determine  the  initial  stressed  state,  it  is  necessary  to  satisfy 
boundary  conditions  (Xa  =  0) 
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0 


7(0 

‘  ijO 


(1.7) 


The  stressed  state  also  depends  on  conditions  as  x^  ->■  ±<*>,  Xa  -*■  ±“.  The 
simplest  stressed  state  results  when  conditions 


limo“,=0,  lim=“,,=  0  (1.8) 

are  satisfied. 

In  this  context,  all  components  of  the  stress  tensor  are  equal  to  zero 
except  Oaa/  which  is  defined  as  follows 


2 


(1.9) 


(for  the  second  and  third  versions  respectively). 

Formulation  of  the  problem  of  wave  propagation  in  the  given  three  ver¬ 
sions  of  the  initial  state  will  depend  on  the  physical  properties  of  the 
material  of  the  half -space .  Here  problems  are  possible  where  the  ma¬ 
terial  is  either  a  dielectric,  or  an  ideal  conductor,  or  a  superconduc¬ 
tor. 


2.  Of  considerable  significance  in  derivation  of  boundary  conditions  is 
allowance  for  the  behavior  of  the  normal  to  the  perturbed  surface  of  the 
half- space  [3]. 


Let  the  perturbed  surface  of  the  half -space  be  defined  by  the  equation 

A',  ■  «,(A',.  0,  X,.  t)  (2.1) 


Then  the  approximate  (linearized)  expression  for  the  normal  to  the 
perturbed  surface  takes  the  form 


n 


dx. 


dx. 


(2.2) 


To  derive  linearized  boundary  conditions  on  surface  Xa 
start  from  the  following  general  conditions;  when  Xa  =  0 


n=TY  ■  n 


0,  We  will 


(2.3) 


where 


3i  =  9ni 


—  A  ^  ^  L 

Ti  -  7  /  li  +  7'j2  /  +  7^»3^ 


(2.4) 


Using  expressions 

for  the  perturbed  state,  and  linearizing  condition  (2.3)  with  allowance 
for  smallness  of  perturbations  and  with  consideration  of  (1. 1-1.6), 
(1.9),  (2.2),  we  get  the  boundary  conditions  on  the  surface  of  the 

half -space. 

In  all  three  versions,,  boundary  conditions  for  components  of  the  per¬ 
turbed  electromagnetic  field  are  the  same,  and  take  the  form 


e,-e\'^ 


2//n  dM, 

r  dt  ’ 


0)  (2.6) 
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Conditions  for  stresses  for  versions  (1.1)-(1.3)  in  the  corresponding 
sequence  are  obtained  in  the  form 


-0,  3„- 

4t: 

h,= 

4- 

.  X. 

333  4-  —  A,  =•  - 
4Tr 

- 

4/^ 

(2.7) 

H\  du. 

Ho 

(2.8) 

2r  dx^  ’ 

!2  —  ^ 

- ?  hit) 

2-  = 

H\du, 

Ho 

HI  dti. 

(2.9) 

2-  dx^ 

■  1 

=0.  + 

~  ^  'dXa 

In  (2.7)-(2.9)  and  below,  the  primes  of  o[.j  have  been  omitted. 

Let  us  take  a  look  at  some  special  problems  of  surface  wave  propagation. 

3.  Let  the  initial  magnetic  field  correspond  to  version  (1.1).  The 
material  of  half-space  Xa  >  0  is  a  dielectric.  The  problems  of  planar 
and  antiplanar  deformation  relative  to  planes  (x^,  Xs)  and  (Xa,  Xg)  are 
separated. 


Let  us  consider  the  plane  problem  where  displacements  and  components  of 
the  perturbed  electromagnetic  field  do  not  depend  on  coordinate  Xa. 


We  give  the  equations  and  boundary  conditions  of  the  problem.  The 
equations  of  the  theory  of  elasticity  in  displacements 


>  'ru. 


o'/  It]  I  ri  ,  .)  .. 

'oxt  ^dx^_  ^  '  -  dx^()Xj 


(rU„ 


2  —i  r  Vi 

'^dx\  ^  ' 


2^2 


dx\ 


dhu 


_ 

d^u„ 


dx^dx^  dt'‘ 


(o.i: 


Equations  of  electrodynamics  in  region  Xa  >  0 

_^Jh.  —  ^  -  _i  ^  ^ 

d.Vj  c  dt  dxi  c  at  ’  dx-^  dx^  c  dt  '  dx^  dx^ 

(3.2) 

in  region  Xa  <  0 


. 1  I  de\p  de<if'> 

dx^  ~  c  dt  '  '■  dJCj  “  c  at  '  d.r,  dx^ 

•  ■  dep:  ^ 

dx^  dx^ 


1  dh<‘^ 
c  dt 

(3.3) 


According  to  (2.6)  and  (2.7),  the  required  boundary  conditions  at  Xa  =  0 
take  the  form 
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.... 

diii 

da, 

aw. 

=0 

c 

"dt ' 

'  dx, 

axj 

(.  -f.  2(/)-f"-i  ..|  / 

dw, 

/lo. 

A*'’^ 

ax, 

ax. 

4t: 

4r 

3 

(3.4) 


Representing  the  sought  functions  in  the  form 


/(Xj,  X;,  i) 

and  following  the  scheme  of  solution  of  the  problem  of  Rayleigh  surface 
waves  [4],  we  get  an  equation  that  determines  the  velocity  of  a  surface 
wave  (it  is  assumed  that  e  =  1,  «  1) 


(2-  r.  i  // 1  ^1(2  -  r,)  -  4/1  -V 1  -  h 


(3.5) 


where 


..  Jl^ 


For  Xi  =  0,  equation  (3.5)  coincides  with  the  Rayleigh  equation.  Just 
as  for  the  Rayleigh  equation,  it  is  shown  that  equation  (3.5)  has  only 
one  real  root  satisfying  condition  0  <  n  <  1- 


Let  us  give  the  equations  of  the  corresponding  antiplanar  problem. 

The  equation  of  the  theory  of  elasticity  in  region  >  0  takes  the  fom 


,2/^.  i 

\dx\  ^  dxl  )  dt* 

The  equations  of  electrodynamics  in  region  Xa  >  0 


(3.6) 


dx^  dXi  c  dt 
in  region  Xa  <  0 


_ ^  ,  1  ^  +  ^-0 

ax,  ~~  c  h'  dx^  t  dt  '  dx,  dx, 

(3.7) 


5x7 


ax. 


1  ^ 

T  ~dT'  ax, 


Boundary  conditions  at  Xa  =  0  take  the  form 


1  1  aA^‘)  dh\'> 

T  W'  lH' 

(3.8) 


A,=Ai'^  O^-f^*A,=  0  (3.9) 

ax,  2it 

We  have  the  task  of  finding  the  solutions  of  equations  (3. 6) -(3. 8)  that 
satisfy  boundary  conditions  (3.9)  and  conditions  of  damping  as  Xg  ±®°. 
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In  solving  the  problem  in  this  formulation,  we  come  to  the  necessity  of 
an  additional  boundary  condition  at  Xa  =  0.  This  condition  is  given  in 
the  form 


1 - 


—  / 
c  OXj 


(3.10) 


which  is  due  to  the  jump  in  the  tangential  component  of  the  magnetic 
field  in  consequence  of  a  change  in  the  direction  of  the  surface  cur¬ 
rent. 


With  allowance  for  additional  condition  ( 3 . 10 ) ,  the  problem  is  complete¬ 
ly  defined.  Solution  of  the  problem  shows  that  there  is  a  purely  shear 
surface  wave  with  velocity  defined  by  the  formula 


(3.11) 


4.  Let  the  initial  magnetic  field  be  assigned  in  form  (1.3)  (third 
version).  It  is  assumed  that  the  ideal  conductor  model  is  applicable  to 
the  material  of  medium  Xa  >  0.  It  is  further  assumed  that  the  sought 
values  do  not  depend  on  coordinate  X3.  In  this  case,  the  problem  of 
plane  deformation  leads  to  solution  of  the  following  eguations  (V2  = 
Ho/irp): 


‘'H 


31/? 


-T-(ci-d  -  1'?) - 

2  ‘  dx^ax^ 


dt* 


(4.1 


dx-^dxj 


'  dl^ 


for  conditions  on  boundary  x^  =  0 


.da,  , 

-’7r,^V^2r.)dx\ 


(4  2) 


The  solution  of  equation  (4.1)  with  condition  (4.2)  and  the  condition  of 
damping  at  infinity  (x^  ->  <»)  shows  that  the  condition  of  existence  of  a 
surface  wave  is  the  same  as  in  the  Rayleigh  problem  (0  <  q  <  1).  The 
equation  that  defines  the  velocity  of  the  surface  wave  takes  the  form 


(2-r  3/2-ri)  {( 1 +6P)  1 2  +  P(-l/2+  40  0,3/2)]  -  ( 1  -0,3  )r,}  = 
=  4(]+9P)3'2(l+303/2)>/2(l+3/2)'/2(l  -3/4){ l--^)''2(  1-0-0'''.  3= 


(4.3) 

V]lc\ 
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Investigation  of  equation  (4.3)  relative  to  the  velocity  of  the  surface 
wave  n  shows  that  the  equation  has  two  real  roots  in  interval  (0,  1). 
One  root  corresponds  to  the  velocity  of  the  ordinary  Rayleigh  wave  with 
allowance  for  the  influence  of  the  magnetic  field.  The  other  root  shows 
the  presence  of  a  surface  wave  with  a  velocity  much  lower  than  that  of 
the  Rayleigh  wave.  This  velocity  can  be  computed  from  the  following 
approximate  formula: 

27:(/-fG) 

The  problem  outlined  in  this  section  has  also  been  considered  in  [5-7]. 
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Diffraction  of  Shock  Wave  Around  Wall  Making  Obtuse  Angle 

917F0082B  IZVESTIYA  AKADEMII  NAUK  AEMYANSKOY  SSR:  SERIYA  MEKHANIKA  in 
Russian  Vol  43  No  1,  1990  (manuscript  received  10  May  88)  pp  10-18 

[Article  by  L.  D.  Azatyan  and  A.  G.  Bagdoyev,  Institute  of  Mechanics, 
ArmSSR  Academy  of  Sciences ] 

[Text]  The  paper  considers  the  problem  of  diffraction  of  a  strong  shock 
wave  around  an  elastic  wall  that  foms  an  oblique  angle.  Lighthill  used 
Chaplygin  transformation  to  study  the  analogous  problem  and  get  an 
analytical  solution  in  the  case  of  a  solid  wall  [ 2 ] . 

§1.  Formulation  of  the  Boundary  Value  Problem 

Let  a  shock  wave  front  move  at  velocity  Uo  through  a  quiescent  gas  and 
impact  at  time  t  0  on  a  corner  with  sides  that  fom  a  small  angle  e 
with  a  break.  Fig.  1  shows  the  pattern  of  motion  for  some  t  0.  The 


Fig.  1 
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parameters  of  homogeneous  flow  (1)  behind  the  incident  shock  wave  are 
determined  from  conditions  of  dynamic  consistency  [1] 


1 


T(T^l) 


( V  •  r  1  )  • ' 

V,--  !)-■!“  : 


f'ln  ‘/i  — 


iM=-: ) 


-,.•1 


^0 


(I.l) 


(designation  of  parameters  is  conventional). 

Let  us  determine  the  flow  in  region  (2)  behind  the  reflected  shock  wave. 
In  the  case  of  an  elastic  wall,  assuming  that  motion  in  the  wall  is 
homogeneous,  we  get  in  the  elastic  region 

d-TV  _  1 

dt'^  ~  >.+2u- 


where  A,  y  are  the  elastic  Lam^  constants,  pcT  is  the  density  of  the 
elastic  body,  w  is  normal  displacement  of  the  elastic  medium,  and  a  is 
the  velocity  of  propagation  of  the  elastic  longitudinal  wave.  On  bound¬ 
ary  X  =  0  from  the  condition  of  concomitant  motion  of  the  wall  and 
nearby  particles  of  fluid  (gas)  we  have 


q.M 


dt 


Y^hgre  n  is  the  normal  to  the  wall ,  and  the  solution  in  the  elastic 
region  is  w  =  f(x  +  at).  Writing  the  conditions  of  continuity 


=?aaV  '(^^0 


we  can  get  two  conditions  on  the  wall 


whence  we  have 


qin=af\at), 


(1.2) 


where  6  =  PzCz/pcrr-a- 

Thus,  accounting  for  a  one-dimensional  model  in  an  elastic  medium  gives 
a  boundary  condition  in  the  form  of  a  Winkler  base.  Flow  parameters  in 
j^egion  (2)  are  determined  from  relations  for  an  obligue  shock  wave  [1]. 
Eliminating  Pz,  Pa  from  these  expressions,  with  allowance  for  (1.2),  we 
find  a  cubic  equation  in  order  e  for  determining  the  normal  velocity  Vo 
of  the  reflected  shock 
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—<•/,)  1/34- lX(aP,-(7i)( I  rg-f  I(-'-3)p^<7J  f- 

-t-  (v  vV-2  ,7'’j(7Pi-v,)  f  - 

-27p=^/t-r2?7o2^’^P,-(Y--l)7^j<^  7=0/p,c',  (1.3) 

Then  in  the  same  approximation  from  expressions  for  an  oblique  shock 
wave  we  get  the  values  of  the  remaining  parameters  for  a  plane  reflected 
shock 

^  ^o)  1  ~f  ^ ())  I 

l+ap,(<7j-t- l^o)  ’  «7^j  -i  «Pi(^/i  l-V'o)*-i-V''o 


<73-.=7,(-’  "®').  ^2 


where  e '  is  the  angle  of  reflection  determined  from  the  condition  of 
inseparability  of  the  incident  and  reflected  shock  waves  at  point  M 

(1.5) 

In  the  region  of  nonuniform  flow  (3),  we  will  give  all  quantities  the 
subscript  "3."  The  boundary  conditions  of  diffraction  are:  arcs  AB  and 
A'B'  that  are  the  fronts  of  perturbations  produced  by  small  angle  e, 
wall  AOA',  and  curved  shock  wave  BCB'  resulting  from  diffraction  away 
from  the  vertex. 

The  problem  is  to  determine  the  parameters  of  motion  in  region  ( 3 ) . 
Since  angle  e  is  small,  quantities 


P.-P* 


P=f'.5-:2-  «(“•  'I')  “93-92 


will  be  of  order  e,  and  the  theory  of  perturbations  can  be  used  to  find 
these  quantities.  Since  there  is  no  characteristic  dimension  of  length, 
the  problem  can  be  considered  self-similar,  and  then  the  equations  of 
gas  dynamics  after  linearization  with  respect  to  e  in  variables 

x^g2xt  y 


are  written  as 


da  da  dP  .dv  ,  dv  dP  .dP  dP  da  .  dv 

+  t 

di  OTi  di  dr^  dri  (?;  di]  <7;  d^ 


On  plane  ri/  the  diffracted  shock  wave  BOB'  and  the  wall  AOA'  can  be 
approximated  by  straight-line  segments  BB'  and  AA'  (Fig.  2). 
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In  this  context,  the  value  of  ^  on  the  wall  is  equal  to  zero,  and  on  the 
curved  shock  the  value  of  C  is 

Arcs  AB  and  A*B*  in  system  n  coincide  with  the  arc  of  a  unit  circle. 

The  solution  of  system  (1-7)  should  satisfy  the  defined  conditions  on 
the  boundary  of  the  diffraction  region.  On  perturbation  fronts  AB  and 
A*B*,  the  functions  are  continuous,  and  consequently  on  these  arcs 

p  =  p=:0,  i;)  =  0  (1-8) 

According  to  (1.2),  we  have  on  the  wall 

^  ~  ^^2  ._9p 

[^2*'  2  I  a  I  *  .  a 

implylnq  that  on  the  wall  in  region  (3) 

c/9-=” 

Then,  as  implied  by  the  second  equation  of  (1.7),  at  C  =  0 

^  _  Or  ^  ^ 

d\  ‘  dr, 

To  find  the  boundary  condition  on  ^  =  ko,  we  represent  the  perturbed 
shock  wave  front  in  the  form  C  =  ko  +  where  i|;(ri)  is  a  function 

that  defines  the  unknown  shape  of  the  shock  wave.  Writing  the  condi¬ 
tions  of  dynamic  consistency  on  the  curved  part  of  shock  BB'  and 
linearizing  them  relative  to  the  flow  behind  the  plane  reflected  wave, 
we  get  the  following  system  of  equations: 


«•  2— (pj-r'l)  l-'o)  =  O'  P'A'^  f  2(>,«(i/2.,—  I'o)  -'Mix—  =  0 

(1.10) 

V.V"  r«lo-  -^C,P-r  -  ‘  ^-72.(«-V")=0 

V  ._  1  V  --  1  0, 


where  the  projections  on  the  normal  and  the  tangent  to  the  shock  wave 
are  written  in  the  linear  approximation  in  terms  of  projections  on  axes 
Ox,  Oy,  Vo  +  V  is  the  velocity  of  the  curved  shock  wave,  V  = 
C2(i|J  - 

The  solution  of  system  (1.10)  takes  the  form 

P=Q(’!>-7ii')/£,  u  =  7iy)/£.  v=:-qi'y-qit^z') 

Q={(p-p2)l2l^2  +  (T-l)(^2.-V^o)M-(T-l)(<7i  +  9^..H92..-V'o)P^)(9-V"o) 

(1.11) 

G={(i-i  )(qx-^qix){qix—VQ)Y2-  T(pi-p»)[^2-^  V'o)*l}<^» 

E  «  p,(  v^f  ] 

Eliminating  function  tp  from  (1.11),  we  get 

U=AnP,  —  on  ;==Ao  (112) 

Ot;  £77] 

where  Aq  =  G/Q,  Bq  qiE/Q. 

Now  by  using  equations  (1.7)  we  can  eliminate  u  and  v  from  equations 
(1.12),  giving  the  condition  on  the  shock  for  P  alone 


dPldl  ■ri(A^-irk,)-8,k,-,i-^ 

dP.dfi  I -HI  ^ 

In  addition  to  condition  (1.13),  a  condition  is  taken  according  to  which 
the  change  in  v  along  the  shock  from  the  center  C  to  vertex  B  is 


fdv  fBodP, 

\—-dr,=  _«  —dr,=q:, 

J  Or:  J  r.  d  n 


B„  dP 


System  (1.7)  can  be  reduced  to  a  single  equation  for  P 


(1.14) 


dl 


dP\  r. 


(1.15) 


The  Chaplygin-Buseman  transfomation  is  used  to  solve  the  formulated 
boundary  value  problem  with  respect  to  the  solution  of  equation  (1.15) 
with  boundary  conditions  (1.8),  (1.9),  (1.13),  (1.14): 
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:  =  rcoS/,  ■»]— rjin'-i,  ~ 


9rj  r. 

Jl _  3=:arctan  — 


(1.16) 


converting  (1.5)  to  a  Laplace  equation  in  variables  p,  p.  Then  confor¬ 
mal  mapping  is  used  [ 2 ] : 


ii. 


=  f  /y,) 


(1.17) 


Fig.  3 

Function  (1.17)  transforms  diffraction  region  AA'B'B  to  a  rectangle 
(Fig.  3) 


1  ,  1  ^  /'O 

- - T‘"7irr 


(1.18) 


Thus,  the  problem  is  to  determine  the  function  P  that  satisfies  the 
Laplace  equation  AP  =  0  in  variables  Xi,  within  rectangle  (12.18), 
and  the  conditions: 


on  AB  ancW'iS'  we  have  P=  —  —0 


.AA'\  =6 cot  y, 

dPidy^ 


(1.19) 


(1.20) 


on  -|^siny,cosyi=^^/4oCos®yx-So  fto-.(l-Ag)''Ml-21) 

In  virtue  of  equation  (1.20),  there  is  a  difference  between  the  given 
problem  and  the  problem  of  [2]  in  which  9=0,  while  the  fairly  complex 

function  a)(Zi)  =  t - iq —  is  determined  by  trial  and  error. 

S2.  Solution  of  Formulated  Boundary  Value  Problem  by  Method  of  Separa¬ 
tion  of  Variables 

We  will  seek  the  solution  of  the  Laplace  equation  AP  =  0  in  the  form  [3] 


AXi,  y,)=7('(Xi)y'()i) 


In  this  case  we  get 
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A‘~.4fxp(v'.A:i)  +  /3exp  (  -•'-’Cl).  )'=Csinv'.y,  l-DcoS/yi 


Satisfying  condition  (1.19),  we  get 

P=v  {Ak^dxp(kx^)  rBk*i'ip(—^x^))ilnky^  (•2.1) 

1 

Condition  (1.20)  can  be  written  as 

V  A|(.4. — Z^*)sin  A’yj  sinyi  — 0( .‘I*--*  /9*)cos  Viccs/cy,]  =0  ('-.2) 

*--^1 

and  condition  (1.21)  can  be  written  as 


S’  k  (.Aie.vpifc' )— -f 


V 


X  (^Aa:<)S'2yi-Ao 


(2.2 1 


Multiplying  (2.2)  and  (2.3)  by  cos  ly,  (1  =  0,  1,  2,...)  and  integrating 
from  -IT  to  IT,  we  get  a  system  of  equations  for  determining  A,^  and  B,*. 
It  can  be  shown  by  simultaneous  solution  of  these  systems  of  equations 
that  Ak  and  B,,  with  even  k  are  equal  to  zero,  while  A,^  and  Bv,  for  odd  k 
will  be  determined  from  the  following  system  of  equations: 


yi,(l_&)=fij(l.+  0).  (l  +  5).i,-(lU)/5,=31(l  -6).43-(1t6)Z?3] 

+  3f>l,(l  +A,)e^'-B,(\-Ao)e-'^']=0  (2.4) 

.4„+3  =  l(l+0)(/^-/?)-(l-A)'Vexp(-(/i+3)).)l/Wf' 

5„+3*I  ( l-6)(/=-->9)-(  l+4o):Vexp((«+3)).)]/  W 


1)1(1— /4o)A„_ie.'cp((«—l)/.)-(l-i-i4u)fl„_  iexp(— (rt— 1)/)] 

/?=2(«+ 1  )^/lo-2  ^)('4„+,exp((«  f  1  )l)4-5„„exp( - («-f  !)>-)] 

A^=(«-fl)I(l+0)'4„+i-(l-e)S„+,l 
«?'  =  (/i+3)[(l  +  ^o)(l+O)exp((/2+3)l)-(l->lo)(l-0)exp(-(«  +  3)/.)l 

From  (2.4),  all  A^^  and  B,^  are  determined  in  terms  of  Aj^  (n  >  2).  Condi¬ 
tion  (1.14)  is  used  to  find  A^.  In  what  follows,  A^,  =  A,^/Ai,  B^:  = 
B,*/Ai.  The  convergence  of  series  (2.1)  is  proved  as  in  [3]. 

S3.  Determination  of  Gas  (Fluid)  Flow  Parameters  in  Diffraction  Region 

We  write  expression  (2.1)  for  perturbed  pressure  in  variables  5,  r|. 

From  formulas  (1.15)  and  (1.17)  we  have 
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Yj  =  arctan 


V 


(3.1) 


■*^1= 


With  the  use  of  (3.1),  expression  (2.1)  in  plane  n  takes  the  form 


p=y 

“o 


X  sin  ('?nX  1  )  -^  r  (2«  t  0 a«±an  ^ 


(3.2) 


Components  u  and  v  can  be  determined  from  the  equations  of  system  (1.7). 
To  find  coefficient  Ai,  we  find  component  v  on  g  =  0  from  the  third 
equation  of  system  (1.7) 


r_L^rfT;=y  (— l)'‘  +  '(2«-f  l)(^.!n^i-rfi--.4-i)  X 
xf  arcsl:i  Tj-i- —  s‘n(2/s  arc  sin  tj)  (3.3) 

'  L  *  -* 

The  quantity  vatC  =  0,  ri  =  0is  equal  to  the  value  of  and  hence 

^  1  S  2/1  + 1 )( -ia, .  I  +  M ) 

/l-O  ^ 

where  is  determined  from  formula  (1.4).  Hence  we  have 


ii 

.•u 

?o 

.M---1.G 

M  -2 

0  --  0 

0-0.1 

f):-0.3 

fi  -vO 

fi=:=0,l 

0=0,3 

fi  0 

0-0.1 

ii=0,3 

P(l) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P(0.5) 

0.24 

0.21 

0.16 

0.41 

0.36 

0.31 

0,49 

0.43 

0.37 

1 

P(0) 

0.28 

0.15 

0.21 

0.-I8 

0.44 

0.40 

p.58 

0.52 

0.47 

H-0.5) 

0.24 

0.21 

0.16 

0.41 

0.36 

0.31 

0.48 

0.43 

0.37 

p(-i) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

PiO) 

Ml 

1.05 

0.94 

1.92 

1.85 

1.74 

2. 68 

2.61 

2.49 

/i(0,5) 

1.07 

l.Ol 

0.91 

1.79 

1,73 

1.65 

2 .50 

2.45 

2.36 

P2(0) 

I  »05 

1.00 

0.90 

1.78 

1.72 

1.64 

2.46 

2.41 

2.32 

Pi  (-0.5) 

1.07 

1.01 

0,91 

1.79 

1.73 

1 .65 

2.50 

2,45 

2.36 

P,(-l) 

l.ll 

l.Oo 

0.94 

1.92 

1.85 

1.74 

2,68 

2.61 

2,49 

15 


_ 2^ _ _ 

^0 


) 


Numerical  calculations  are  given  to  clarify  the  qualitative  picture  of 
reflection  of  a  shock  wave  from  the  wall.  From  formula  (3.2),  the 
pressure  on  the  wall  is  computed  at  Mach  numbers  M  =  Vo/cq  =  1-2,  1.6,  2 
and  0  =  0,  0.1,  0.3.  Calculations  for  the  distribution  of  dimensionless 
pressure  P  =  -P/eUo  on  the  wall  show  that  the  pattern  that  occurs  for  a 
rigid  wall  is  qualitatively  retained,  and  moreover  there  is  a  decline  in 
pressure  P  with  decreasing  6.  The  distribution  of  total  pressure  P3  on 
the  wall  is  also  calculated  for  e  =  0.1,  pressure  P3  decreasing  with 
increasing  pressure  0.  The  table  also  shows  that  P  and  P3  on  the  wall 
in  the  diffraction  region  increase  with  an  increase  in  the  Mach  number 
at  fixed  0. 


M={,2  MM.6 


Fig.  4 


Fig.  4  shows  graphs  of  the  distribution  of  pressure  P  on  the  wall  for 
two  cases. 
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[Article  by  V.  Ts.  Gnuni,  Institute  of  Mechanics,  ArmSSR  Academy  of 
Sciences] 

[Text]  In  the  linear  theory  of  bending  of  plates  under  a  transverse 
load,  it  is  known  that  forces  (1/  k  -  1,  2)  do  not  arise  in  the 

middle  plane  of  the  plate,  and  hence  it  is  assxamed  that  the  transverse 
load  acting  on  the  plate  does  not  affect  frequencies  of  free  oscilla¬ 
tions.  However,  in  bending  of  very  thin  (flexible)  plates,  considerable 
tensile  and  shear  forces  arise  in  the  middle  plane  of  the  plate  for 
flexures  comparable  with  the  thickness  of  the  plate,  and  these  forces 
may  significantly  alter  the  frequencies  of  free  oscillations  of  the 
plate . 

This  paper  attempts  to  give  a  qualitative  picture  of  the  influence  of 
transverse  loading  on  the  frequencies  of  free  oscillations  of  a  plate. 

I.  Let  a  flexible  rectangular  plate  of  dimensions  a,  b,  h  be  hinged 
around  the  edge  and  subjected  to  the  action  of  a  uniform  transverse  load 
of  intensity  g. 

The  equations  of  bending  of  a  flexible  plate  referenced  to  a  rectangular 
cartesian  coordinate  system  Oxyz  are  [1] 


— 

Eh 


dx^ 


IF 


dxdy  dxdy 


dx~ 


(l.l) 


where  Wq,  $o  are  functions  of  flexure  and  forces  of  the  initial  static 
state,  D  =  EH^/12(1  -  v^),  E  is  Young's  modulus,  v  is  Poisson's  ratio. 
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The  forces  in  the  middle  plane  are  defined  as  follows  in  terms  of  func¬ 
tion  $o5 


JO  75=-^^,  7';,=  --^  ■  (12) 

”  dy-  dx^  '■  dx  dy 

Small  transverse  oscillations  of  the  plate  about  the  bent  equilibrium 
state  (1.1)  are  described  by  the  equations  [2] 


Eh  dy*  dx* 


d*w  _ 
dE 


_yo 


+  - 

d*w 


-27?. 


dv*  dx* 


d*w 

^  d*WQ  d*w 

^  ’ 

dxdy  dxdy 

,  dht' 

d*w 

^  dx  dy 

d*'\' 

d*<\> 

dx  dy 

dx^  dy^ 

(1.3) 


where  w,  $  are  the  functions  of  flexure  and  forces  of  the  perturbed 
(oscillatory)  state,  p  is  density  of  the  plate  material,  and  t  is  time. 

Equations  (1.3)  account  for  both  the  forces  and  flexures  Wq  of  the 

initial  state. 


Equations 


ch,,=:ct)(,„slii  A, .x: Sill  I4y.  sin  s-vsin  UiV 

f'.  sill  s!n  !'.„y  cos  <'.„,„l!,  7C>— ;',m  sin  ■■  r-.X  sin  JinV  cos 


(1.4) 

(1.5) 


where  =  mir/a,  =  nir/b  identically  satisfy  the  boundary  conditions 
for  hinged  support  of  a  plate 


“ll 

p 

ii 

O 

ic'^O, 

.14,, =0 

for 

11 

p 

11 

(1.6) 

^32—0, 

u—O, 

4422=0 

for 

y=0,  y=:* 

As  shown  in 

[1] 

,  expression 

(1.4) 

within 

the 

limits  of  flexures  of 

initial  state 

Wo 

of  the 

order 

of  two 

or  three  plate  thicknesses 

gives  a 

rather  good  description  of  the  initial  bent  state. 

For  (1.4),  (1.5),  using  the  Bubnov-Galerkin  method,  systems  of  equations 
(1.1),  (1.3)  yield 


(1.7) 


(1.8) 


^nfon^^uJ^ 


16 


on  2 


_  /i  , 

nn  —  I  / 
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Here  we  have  introduced  the  notation  to  denote  the  frequencies  of 

free  oscillations  of  the  plate  about  the  bent  equilibrium  position,  and 


(1.9) 


to  denote  the  frequencies  of  free  oscillations  of  the  loaded  plate 

Eh  '  ''  a  1 


Knin  —  tn  '  ^mn  — 


32^2 


m’.v. 


3 


("*L 

V"l'm  / 


(I.IO) 


1 


a^h‘‘  (4.'7<“  — 1 '  (4«"— 1 ) 

As  can  be  seen  from  (1.8),  the  frequencies  of  free  oscillations  of  a 
plate  loaded  by  pressure  q  increase  as  compared  with  a  plate  free  of 
pressure,  and  clearly  the  influence  of  the  initial  state  can  be  disre¬ 
garded  for  flexures  of  the  initial  state  that  are  much  less  than  the 
thickness . 


2.  By  way  of  illustration,  consider  the  example  of  a  plate  of  square 
planform  at  v  =  0.3.  In  this  case,  equation  (1.7)  is  written  as 

yf+2.51/*-v»=0  (2  1) 

where  f*  =  fon/h,  q*  =  9a'‘q/8tr^h'^E  are  dimensionless  parameters  of 
deflection  and  loading. 

TABLE  1 


1 

1 

2 

i 

1 

3 

4 

5 

10 

15 

20 

I  * 

0.37 

0.6/ 

i  *1 

\f2 

1.8 

2*1 

2.4 

Table  1  gives  values  of  f*  for  several  values  of  q^. 

From  (1-8)  for  frequencies  of  free  oscillations  of  a  loaded  square  plate 
at  V  =  0.3  we  have 


.  ''  l/- 


~  |/  ~  -  rrf  'r0,2( 


(-:.2) 


where  we  have  introduced  the  notation 


_  m-nHm“+n^-0^5) 

(,^0.25)(«^-0,^) 

On  the  basis  of  (2.2)  for  various  values  of  q^  we  determine  the  normal¬ 
ized  frequencies  of  free  oscillations  of  a  loaded  plate; 

rt"  I  /  T  (2.3) 
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for  m  =  n  =  1  for  which  (!)„„  - >  min. 

cn  »  n 

TABLE  2 


0 

1 

2 

3 

4 

5 

10 

15  ' 

20 

*"11 

2 

2.6 

3.3 

3.9 

4.6 

5/2 

7.1 

8.4 

9.3 

Clearly,  accounting  for  forces  and  flexures  of  the  initial  state  has  a 
fundamental  effect  on  the  frequencies  of  free  oscillations  of  the  loaded 
plate  for  flexures  comparable  with  thickness. 

The  numerical  results  obtained  in  section  2  are  general  for  arbitrary 
plates  of  square  planform.  However,  it  can  easily  be  shown  that  these 
results  are  competent  only  for  sufficiently  thin  plates  that  allow 
bending  comparable  with  the  thickness  of  the  plate  within  the  limits  of 
elasticity  or  up  to  brittle  fracture.  From  the  condition 


.2 


where  fundamental  stresses  o 


l  -t  i  —2 

n  ®11"22  »  ^2 

2  + 

(2.4) 

^  f  ^22?  i 

r  ®  .12 

are  defined  by  formulas  [1] 

_L  <^"^0 

Ez 

/<)=Wo 

,  ,,  ^^0  \ 

h  dy* 

1-/ 

V 

^  df  ) 

/d-Wa  ^ 

\ 

1  T 

I 

Ez 

d^w 

\  dy^ 

dx^  ) 

1.  ^12 

h  dxdy 

1  +  v 

dxdy 

we  can  get  the  following  estimate  for  the  admissible  relative  thickness 
h/a: 


h 


Here  f^  for  prescribed  q  is  determined  from  equation  (2.1). 

Thus,  for  prescribed  a,  E,  there  are  real  plates  for  which  ac¬ 
counting  for  the  initial  state  is  significant  in  determining  the  fre¬ 
quencies  of  free  oscillations.  For  exaii^Jle,  at  a^^/E  =  0.005  and  f^  = 
2.1,  a  =  200  cm,  we  get  h  i  2.53  cm.  Let  us  note  that  f^  =  2.1  corre¬ 
sponds  to  a  normalized  load  of  q^  =  15. 

In  conclusion,  we  note  that  questions  of  the  influence  of  an  external 
statically  applied  load  on  the  frequencies  of  free  oscillations  have 
been  rather  thoroughly  studied.  And  in  the  case  of  shells,  the  initial 
state  caused  by  external  loads  may  be  treated  as  zero-torque  [3]  or 
significantly  torqued,  as  well  as  nonlinear  [2,  4].  On  the  other  hand, 
in  the  case  of  plates,  only  geometric  nonlinearity  enables  accounting 
for  the  effect  of  the  initial  state  on  frequencies  of  small  oscillations 
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about  the  bent  equilibrium  state.  In  this  sense,  we  can  derive  rela¬ 
tions  for  studying  the  frequencies  of  free  oscillations  of  plates  of 
rectangular  planform  and  for  circular  plates  loaded  by  transverse  pres¬ 
sure  from  the  results  of  [2,  4]  obtained  for  a  cylindrical  panel  [2]  and 
for  corrugated  shells  of  revolution  [ 4 ] . 

Let  us  also  note  that  we  can  derive  an  approximate  analytical  expression 
on  the  basis  of  [2]  for  the  frequency  of  free  oscillations  of  a  loaded 
plate,  and  [4]  can  be  a  basis  for  exact  numerical  analysis  of  frequen¬ 
cies  of  a  circular  pressure-loaded  plate. 
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[Text]  Consider  the  problem  of  designing  an  orthotropic  circular  plate 
of  radius  R  and  thickness  h  under  constant  normal  load  q. 

The  plate  is  referenced  to  a  cylindrical  coordinate  system  0^  r,  0,  z  so 
that  coordinate  plane  z  =  0  coincides  with  the  middle  plane  of  the 
plate.  The  material  of  the  plate  has  three  planes  of  elastic  symmetry 
that  are  perpendicular  at  each  point  to  the  corresponding  coordinate 
lines,  the  pole  of  anisotropy  coinciding  with  the  center  of  the  plate. 
The  plate  of  greatest  stiffness  is  determined  by  varying  the  structure 
of  the  material  and  the  location  of  the  support  contour. 

It  is  assumed  that  the  plate  is  made  up  of  2n  layers  alternately  rein¬ 
forced  in  the  radial  and  circular  directions.  In  this  arrangement,  the 
layers  with  odd  numbering  (1,  3,...,  2n-l)  are  made  up  of  m^.  elementary 
layers  that  are  reinforced  in  the  radial  direction,  while  the  layers 
with  even  numbering  (2^  4,...,  2n)  are  made  up  of  m^  elementary  layers 
reinforced  in  the  circular  direction.  From  the  foregoing  we  see  that 
the  total  thickness  is  h  =  nm6o/  where  6o  is  the  thickness  of  a  mono- 
layer,  and  m  =  mi  +  m^.  The  parameter  g  =  m^/m  defines  the  relative 
thickness  of  the  radially  reinforced  layers  of  the  stacked  plate.  It  is 
assumed  that  elementary  layers  reinforced  in  both  the  radial  and  circu¬ 
lar  directions  have  the  same  content  of  reinforcing  material  per  unit  of 
volume,  so  that  the  elastic  characteristics  of  these  layers  in  orthogo¬ 
nal  directions  are  the  same  [ 1 ] . 

Let  the  plate  be  supported  around  line  r  =  Ri,  where  0  <  Ri  ^  R.  The 
equation  of  equilibrium  of  an  orthotropic  plate  relative  to  bending  W(r) 
is  represented  as  [2] 


22 


(1) 


dnv<  1  J. 

Drr-^-2D,r^  r''  c/r*  /■’ 

where  i  =  1,  2  for  regions  0  S  r  ^  R  and  R,  ^  r  S  R  respectively. 
Boundary  conditions  are  written  as 

llm  U",(/-)<oc,  0  for  r -A) 

r^o  dr 

ir,=  u;=u.  — ^  =  for  r=R, 

dr  dr 

/V/(-'>  =  0,  ,V<->— (')  for  r—R 


"  dr'- 


,,  /dMV.  I 

('i==  —OrrX - - -  riJl-i—  —r 

V  d,-*  r  dr-  )  /"*  «r 


t  J'.r, 


are  the  bending  moment  and  the  transverse  force  in  cross  sections 
r  =  const  respectively,  and  D^,.,  Dee  and  D^-e  are  the  stiffnesses  of  the 
plate  of  given  structure  [1] 

\^B B,.A\ -  B.-.V  _\KB.,-B^,\,--B.,\n-:  ^ 

Drr-^- - - - -  - - r;  ■  ,.T 

12  *- 

Here  Bee  and  B,^e  are  the  elastic  characteristics  of  a  monolayer  of 

orthotropic  composite  material. 

In  deriving  these  expressions,  terms  of  order  1/n  and  1/n^  were  ignored, 
where  n  is  the  total  number  of  monolayers  reinforced  in  the  radial  and 
annular  directions. 

The  solution  of  boundary  value  problem  ( 1 )  and  ( 2 )  with  allowance  for 
(3)  takes  the  form 


Tf/ _ dR*  •  t2  _ o2 _ j  ^ 

(l-Brr)i  +  Brr  1  ‘ 


34-a 


qR*  -  12  32 

'64Brrh^  (l-BrrT^  +  Brr 


4(9  -/c“) ! 


Ot  •  2 

2(l-k*)  ' 


I  3+h 

I  {9-k’‘){l-- k){k  i  [A  ~ 


2k(\-k^)ii  rk){k  ^n)  I 
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where  the  notation 


ife-‘  •=  %  (  B„i-\  );  r  1 

^rr  (  1  —  -1-  BhU 


-r/R,  %=R,!R 


D,r  (  1 —  Bu:)\  -r  Bh\i 

has  been  introduced,  in  which  B,,e  =  and  Bee  =  Bee/B^,.. 

Our  problem  is  to  find  the  optimum  value  of  parameters  ^  and  p  that 
maximize  the  stiffness  of  the  plate,  i.e. 


with  constraints 

(/-I);  (/=2)  (8) 

and  for  assigned  R,  h,  B„,  Bee,  B^e  and  q. 

By  way  of  exait5)le,  consider  a  plate  made  up  of  monolayers  of  orthotropic 
boroplastic  reinforced  in  the  radial  and  annular  directions. 

An  algorithm  was  developed  for  numerical  solution  of  problem  (7)  with 
constraints  (8),  and  the  corresponding  program  was  written. 


j 

1 

1 

w 

0.00 

0*693 

1*00 

0*6452 

0*05 

0.691 

1.00 

0.5463 

0.10 

0*689 

1.00 

0*4834 

0'15 

0*688 

0*00 

0.4325 

0.20 

0.687 

0»00 

0.3896 

0*25 

0.685 

1*00 

0,3688 

0*30 

0*684 

1*00 

0*3437 

0*35 

0.b83 

1*00 

0*3235 

0*40 

0.682 

1*(0 

0*3072 

0*45 

0,681 

1  I  *00 

0*294 

0*55 

0.680 

0.00 

*•.2674 

0*60 

0.679 

0*00 

0  *2550 

0*65 

0.678 

1.00 

0.J486 

0  *7f 

0.677 

1*00 

0*2442 

0*75 

0,676 

1.00 

0*2410 

0*80 

0.676 

0*00 

0*2332 

0*85 

0.675 

0*00 

0*2262 

Of90 

0.674 

0*00 

0*2211 

0*95 

0.673 

0*00 

0*2195 

l*00 

0,672 

o*co 

0*2262 

Values  of  -  P*,  the  value  of  y  =  _r/R  where  max  Wjl(y,  P,  5)  is 

reached,  and  the  corresponding  min  max  Wi(Y,  P,  g)  are  given  in  the 

P  TT 
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table  for  several  values  of  structure  parameter  g  characterizing  the 
anisotropy  of  the  plate.  Here 


U7=\V'.64/?r./i7l2v^'‘ 

is  is  the  dimensionless  flexure  of  the  plate. 

As  we  see  from  the  table,  the  optimum  value  of  parameter  p  that  charac¬ 
terizes  the  position  of  the  supporting  edge  for  all  values  of  structure 
parameter  C  varies  monotonically  in  the  narrow  range 

0,672  <3  0.693 

The  optimum  plate  design  is  realized  for 

;=0,95,  ^=0,673  and  here  \l''Qp^f=0,219S 

this  design  containing  5%  of  monolayers  reinforced  in  annular  direc¬ 
tions.  The  worst  design  is  realized  for  4  =  0  (plate  made  up  of 
monolayers  reinforced  only  in  the  annular  direction) ,  and  here 

max  min  max  H7=:).6452 

Thus,  structural  optimization  leads  to  considerable  improvement  of  the 
design  of  a  plate  in  the  sense  of  stiffness, 

imix  mill  max  v4^=2,9  i  min  min  max  vi 
i.e.  roughly  by  a  factor  of  3. 
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[Text]  The  paper  examines  the  axisymmetric  problem  of  elastic  wave 
propagation  in  a  radially  inhomogeneous  elastic  ball  with  the  assun^Jtion 
that  radial  and  tangential  stresses  on  its  surface  are  known.  Standing 
oscillations  of  a  radially  inhomogeneous  elastic  ball  are  also  studied. 

1.  General  formulas.  Let  the  Lame  coefficients  A,  p  and  density  p  of 
an  elastic  ball  be  assigned  in  spherical  coordinates  (r,  6,  qj)  by  func¬ 
tions 


4 


where  Aq  -  (y  “  2)Po;  Po»  Po  are  the  Lame  coefficients  and  density  of 
the  medium  at  r  =  ro;  y  =  (Ao  +  2po)/Po. 


The  wave  field  of  perturbations  in  an  inhomogeneous  elastic  ball  r  <  ro 
at  any  time  t  >  0  is  determined  by  some  equation  of  elasticity 


(A-|-2p)grad  div  «— p  curl  curl  «-4-gracl  /  div  «  t 


-  -  d^u 

-t  grad|iX  curl  u  +  2grad  ^  grad  u=o  ~ 


(1.2) 


under  conditions 


«|<=o 


da 

dt 


=  0 


(1.3) 


4-0 


^rr(fu,  0»  8.  ^) — t)\  'r«(r’oi  8,  ^) — 0  (1.4) 


26 


Assuming  Rayleigh's  hypothesis,  according  to  which  the  relation  between 
components  of  the  stress  tensor  and  displacement  vector  Ue,  u^)  of 
an  inhomogeneous  elastic  medium  is  the  same  as  for  a  homogeneous  medium, 
we  have  in  the  axisymmetric  case 


where 


1 

rsir.O 


—  (sin  ^U!,  ) 


(1.8) 


The  displacement  vector  d  of  the  radially  inhomogeneous  medium  is  repre¬ 
sented  in  terms  of  generalized  potentials  as  follows: 


_ 

u(r,  0,  /I -grad  evO+(— )  •“"cur*  curl(0,fr.  9.  0?rl-+-‘^urltl>,(r.e,  t)e\\ 


(1.9) 

i.e.  the  components  of  the  displacement  vector  3  are  determined  from  the 
formulas 


1  d 

= 

sin9  d9 

d9)| 

1 

^  .  "t 

d^dj; 

o 

11 

R  dR 

dRo<i 

ro 

R  dS 

(1.10) 

(1.11) 

(1.12) 


where  is  the  basis  vector  with  respect  to  the  r-axis,  R  -  r/ro,  $*  - 
i  =  1,  2,  3. 


Then,  representation  (1.9)  of  the  displacement  vector  permits  separation 
of  the  vector  equation  of  motion  (1.2)  into  the  following  independent 
scalar  equations  [ 1 ] : 


4  dO; 

\ 

1  d*fl),‘ 

—  0 

(•i-2)/?  dR 

(7-2)7?^  * 

7/?>  d.* 

2v  dch- 

1  rj*(p: 

—  0 

(y-2)/^  dR 

(7-2)/?''  - 

R^  d-* 

k”l 

'  (j)* 

(Y-2)/?  dR 

R^  dt‘ 

(1.13) 

(1.14) 

(1.15) 
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where 


'R'dR'^^dR^  /?••“  sin  0  ^  W/' 


2.  Wave  propagation  in  an  inhomogeneous  elastic  ball  when  its  surface 
is  acted  on  by  radial  time -dependent  pressures.  In  this  case,  gen¬ 
eralized  potential  depends  on  R  and  x,  $^  =  $*  =  0.  Wave  propagation 
in  the  ball  is  described  by  the  equation 


:>7 

under  conditions 


(7 


I 


-  D2 

1 


-0 


dfh, 


iIt-O 


=.0 


^0 


4- do: 


7/?  d,v  I 


■4fp) 

I'o 


(2.1) 


(2.2) 


(2.3) 


Equation  (2.1)  and  condition  (2.3)  in  the  Laplace  transform  with  respect 
to  T 


are 

/7)=  1  -)exp(-/7T)r/- 

as  follows :  ^ 

(2.4) 

,  2y  /p^  .  4  \7r: 

dR-  •  (■  -2)R  dR  \ -j  '  ■;-2//2* 

(2.5) 

1  dR'  dR  l/r.,1 

(2.6) 

The  solution  of  Euler's  equation  (2.5)  under  condition 
boundedness  in  the  center  of  the  ball,  takes  the  form 

(2.6),  assuming 

where 


^I(/?.  P)^  - 


=of{p)R^-^ 

7!^o(’— 


(2.7) 


The  radical  branch  a  is  fixed  by  the  condition  arg  a  =  0  for  p  >  0. 
According  to  ( 1. 10)-( 1. 12)  and  (1.5)-(1.8),  The  components  of  the  dis¬ 
placement  vector  and  stress  tensor  take  the  form 
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where 


Ur{R,  P)=- 


v;i„(a— rf) 


■3v,[R>p)  ■  ^x.AR,  p)——^-~'^of{p)(^^^  " 


Co= 


6--; 
2(t-2)  ’ 


b 


2(3-.-4) 

■.(7-2) 


We  get  the  inverse  Fourier  transform  of  images  (2.8)-(2.10)  under 
tion  f(p)  =  1,  i.e.  f(j:)  =  6(t),  where  6  is  the  Dirac  function, 
general  case  for  any  f(p),  we  find  the  inversion  of  expressions 
(2.10)  from  a  convolution  formula. 

From  formulas  [2} 


a/ -  oA  =  A')  j* 

0 


we  get 


I 


(2.8) 

(2.9) 

(2.10) 


condi" 
In  the 
(2.8)- 


(2.11) 


(2.12) 

(2.13) 


29 


where  H  is  the  Heaviside  function,  is  a  Bessel  function. 

Waves  of  displacements  and  stresses  propagate  in  a  radially  inhomo¬ 
geneous  elastic  ball  when  its  surface  is  acted  on  by  pressures  that 
depend  only  on  time.  The  front  of  these  waves  at  any  t  >  0  is  defined 
by  the  equation 


Equation  (2.15)  defines  a  sphere  with  radius  R  =  exp(-i/YT).  As  a  result 
of  rapid  decay  of  the  radius  of  the  wave  front  in  time,  the  displacement 
and  stress  waves  have  an  appreciable  effect  near  the  surface  of  the 
ball.  These  waves  do  not  reach  the  center  of  the  ball,  and  consequently 
no  reflected  waves  originate  from  the  center. 

In  formulas  (2. 12) -(2. 14) ,  the  terms  associated  with  the  Heaviside 
function  describe  waves  that  do  not  originate  in  a  homogeneous  medium. 
It  propagates  with  a  wake,  and  its  occurrence  can  be  attributed  to 
dispersion  of  transmitted  waves  in  the  inhomogeneous  medium.  Near  front 
(2.15),  the  so-called  wave  with  a  wake  has  a  finite  discontinuity  im¬ 
plied  by  formulas  (2.12)-(2.14)  if  it  be  considered  that  the  limit 
z~^J3(bz)  as  z  ->■  0  is  finite. 

Let  us  consider  steady-state  waves  in  a  radially  inhomogeneous  ball. 
For  boundary  condition  f(T)  =  0oe^‘“'‘  we  get 


Ur{R,  -)=- 


where  phase  U  is  defined  by  the  relation 


(2.16) 


The  phase  velocity  is  determined  from  the  condition  =  const: 


(2.17) 
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(2.18) 


d{\—R) 


The  wave  number  is 


y.  R 


1  1  .  I')*  o  i 


For  the  group  velocity  we  get  the  expression 


<iio  R'l  1  .  ID*  «  I 

dK  *  I  1  “ 


(2.19) 


(2.20) 


As  we  see  from  formulas  (2.18)  and  (2.20),  the  phase  and  group  veloci¬ 
ties  depend  on  frequency  u),  coordinate  R,  and  the  quantity  y. 

Standing  waves  appear  in  a  radially  inhomogeneous  ball  at  frequencies 
lower  than  the  limit 


(2.21) 


3.  Resonant  frequencies  of  a  radially  inhomogeneous  elastic  ball  in  the 
case  of  standing  waves.  Consider  the  solution  of  system  of  equations 
( 1. 13)-( 1. 15) ,  assuming  that  the  generalized  potentials  take  the  form 


cl>*(/^,  'j.  -)=-?i(R)Pnicosij)e-  -,  j=l,  2.  3 


where  P„  cos  6  is  a  Legendre  function,  and  radial  and  tangential  stress¬ 
es  on  the  surface  of  the  ball  are  equal  to  zero. 

Substituting  (3.1)  in  (1.13)-(1.15)  and  (1.4),  we  have  system  of  equa¬ 
tions 

d/T  ^  '2)  R  dR  [  ^  -  *  * 


1 ^  -  I  n(n  \- 1 )  r  i  S  ^ 

^  -  I  mn-fl  rr  =  ^ 


with  boundary  condition 


1  rJh\  1  dllr  .  dUl,  .  =0 


dR  R  11«=> 
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where 


iir^rQ 


^^-rn{n+\)k 
dR 


i  dP„(cos(i) 


Uh — r 0  j  _  ^ 

k 


dh 


/'n 


— -<r. 


rfPn(OOsG) 

de 


^3.8) 

(3.9) 

(3.10) 


The  solution  of  Euler  equations  (3. 2) -(3. 4)  under  condition  of  bounded¬ 
ness  in  the  center  of  the  ball  is  as  follows: 

'?,(/?)  =  A„R^--^;  <f,iR)=Br.R‘  '‘  ;  UR)  (3-  H ) 

where  A„,  B^,  C„  are  arbitrary  constants. 


7. 


Cl: 


Radial  oscillations  of  an  inhomogeneous  elastic  ball  occur  at  n  =  0. 
Solution  (3.11)  at  n  =  0  with  boundary  condition  (3. 5) -(3. 7)  gives  the 
equation  of  frequencies 


//  ( /7  ]  ) 


/I  ' 


T  1) 


The  roots  of  equation  (3.12)  that  define  the  resonant  frequencies  of 
radial  oscillations  are  as  follows: 


v'fl~^  )(’(i— rf)  -  0 


(i.l2) 


Roots  of  frequency  equation  A„(<i))  =  0 


n 

0 

1 

2 

O 

4 

5 

6 

2»5 

2 .3665 
4,4722 

0,0016 

2,8647 

4,4722 

1 ,2245 

4  0.385 
4,4722 

1»6942 

4  r  4722 

5 » 4739 

2,2565 

4,4722 

6,9535 

2,9270 

4.4722 

8  ,.3666 

3.6715 

4,4722 

3 

2,5821 

.3,4642 

2,5821 

3,4642 

0.0030 

3.0817 

3,4642 

1 ,3040 
3,4042 
4,1854 

1,0568 

3,^1642 

5 >477 7 

2,6994 

3,4642 

3.4642 

3,5126 

10 

2,2362 

3,2849 

0,0(1.33 

2,2362 

2,6439 

1.4912 

2,2362 

2.2362 

2.4089 

2,2362 

3,3412 

2,2362 

4,2849 

2 ,2.362 
5,2334 
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(;>.i3) 


At  n  =  0,  substituting  (3.11)  in  boundary  condition  (3.5)-(3.7).  we  get 
a  homogeneous  system  of  linear  algebraic  equations  relative  to  A„,  B^, 
C„.  The  condition  under  which  the  homogeneous  system  has  a  nontrivial 
solution  is  as  follows: 

A„(,n)^{(, -r)I  -  (a, -c)  -  - 2)«(/H  I  )}l(3.  ^  c)0-c-2)  -f 

-  1)(3<„-c--1)(3„— C- 3)  =0  (3.14) 

The  roots  of  equation  (3.14)  define  the  dimensionless  resonant  frequen¬ 
cies  of  a  radially  inhomogeneous  elastic  ball.  They  are  given  in  the 
table  for  several  values  of  n  and  y. 
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[Article  by  Zh.  G.  Apikyan,  Institute  of  Mechanics,  ArmSSR  Academy  of 
Sciences] 

[Text]  In  [1-3],  researchers  consider  introduction  of  a  cylindrical 
solid  into  an  ideally  plastic  tube  ( isotropic  and  orthotropic  cases ) . 

Our  paper  takes  a  look  at  introduction  ( internal  and  external )  of  a  near 
cylindrical  solid  into  a  cylinder  made  of  material  with  power-law  hard¬ 
ening. 

I-  Let  a  cylindrical  tube  made  of  material  taken  as  incompressible  and 
power-law  hardenable  fit  tightly  in  an  absolutely  rigid  cylindrical 
®old.  The  inside  and  outside  radii  of  the  tube  are  a  and  b  respective¬ 
ly.  A  cylindrical  tube  of  considerably  harder  material  with  varying 
outside  radius  R  =  R(z)  =  a  +  vu^  exp  (vz/b),  where  v  and  Ui  are  pre¬ 
scribed  positive  constants  (Fig.  1),  is  coaxially  forced  into  this  tube 
(internal  introduction).  The  material  of  the  second  tube  is  taken  as 
nondef ormable . 

The  problem  is  solved  in  the  cylindrical  coordinate  system  r,  0,  z  fixed 
to  the  rigid  tube  so  that  plane  z  =  0  passes  through  the  entrance  end 
cross  section,  and  the  positive  direction  of  the  Oz  axis  is  opposite  to 
the  direction  of  motion  along  the  axis  of  the  tube.  The  end  z  =  1  of 
the  tube  is  assumed  to  be  free  of  external  forces.  Then  the  components 
of  displacement  u,  w  in  directions  r,  z  do  not  depend  on  0,  and  v  =  0. 

The  equations  of  equilibrium,  relations  between  the  components  of 
strains,  displacements  and  stresses,  and  the  law  of  hardening  take  the 
form 
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where 


^  3ei). 


‘  i  v;,  are  intensities  of  stresses  and  strains, 
I  ,  , 

-  - (~-r  <■  3;,  r  J-) 


We  seek  the  components  of  displacements  in  the  form 


It  ---=  t.f 

where  f  =  f(r),  ^  =  (rf)',  A 
derivative  with  respect  to  r. 

Then 

'Ir 


=  •-/' '  exp('-). 

V 


.‘xp(/.r),  K'  -^e.';p('  :)  ' ' 


is  a  constant,  and  the  prime  denotes  the 


s  i,---  ~evp(  r) 


(1.5) 


0!  -  I 


L=-li^ea'  fxp((w— 1)'~).  ^  '-xpt.'.r) 

I  > 
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where 


m  +  \ 

2'"  3  '  k 


From  relations  (1.2)  and  (1.5)  we  get 

30  =  3,— /.AjU)  — ^expCm'/.r),  3- = — /.tivo  {t  ' (m!  ~) 

x„  =  0.5A',iurexp(m/.-:)  (1-6) 


where  w  =  a)(r)  =  eo"" 

From  (1.1)  and  (1.2)  we  get 

r 

3,  =//(2')-/*jexp(/n/-~)  juj  [  (1.7) 

a 

From  relations  (1.1),  (1.6),  (1.7)  we  find  the  unknown  function  H(z), 

and  for  function  f  we  get  the  fourth  order  differential  equation 


7/(z)=>  .feiCi  exp  (m/  -r) 


(«.  T )"  -  mi*m  (^m.T  +  2  “  Sm*^***  +  (~)  = 


where  is  a  constant. 

Boundary  conditions  take  the  form 

Z)  •,  '^)  (1.9) 

u{b,z)=Q.  «(a,  2:)=v«jexp(v.r  Z»)  (1.10) 

where  and  Pa  are  coefficients  of  friction.  Moreover,  the  condition 


I  3j(r,  l)rdr  —  ()  (1-1 1) 

a 

must  be  satisfied  on  end  z  =  I. 

We  introduce  dimensionless  variables 

a,  .  ,  r  a  .  z  .  I 

Uy  •!  --  l-b  ,  0—-—,  _r>y  —  .  ,  — ,  — 

b  b  b  b  b 

u(r,  z)  ==  ',),  «)(r,  z)  -=  bw*([.,  1),  -rz{r,  -)  { 1. 12) 

3r(r,  2)=fe,a;(p, ',),  ce(/-,  c)  =  a:i3J(p,  ',),  ^^(r,  z)=k^:!]{o,  \) 
fir)  =  by  R{z)^bR*{\).  C,==<bq 
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Fourth  order  differential  equation  (1.8)  can  be  replaced  by  a  system  of 
four  second  order  differential  equations.  In  dimensionless  quantities, 
it  takes  the  form 


where 


df*  _ d'i*  _.l-r4''V  _07» 


uo  2,0 

dF 


do 


2.0 


3 

==  rnFui*  l  ml  -  - - - 

do  \ 


(1.13) 


1  /  /T**  3f  ■i-o2tp*'\tll 


(1  -m)v=p(2vy-(-3)/*?*r'* 


dr*  _  1 

dp  “gn  ‘  \  4  ■  4p 

-  3(5  -3/n)vY**r*— (3-OT)vyf  ♦»r*+4(2/7z— 1  y/Y-p*  r*'‘- 


—  2p*7'*^  -i-4/K'j^p'f*(3/***  7-p*<fi**)  ^ 

—  =  9 /*  *'  -f- _  T'  *  T'  -  v-y  * _ /  ^ -  '  j 

o’  Vo/ 


0 


in~\ 


*^=«)  ^  e*,  =  e,=: 


?jr*‘  +  'FY*  T* 
4o'-  ^  4 


(1.14) 


F*  =zbF  =  -2///v=<f*o.*  -I-  7»«)*/p 

Q*  =  /,'  -’g  =  2pi  r-{  -,f  '  p'‘f  *2)  1  inp-  /  *=) 


Using  relations  (1.6),  (1.7),  (1.12),  (1.14),  we  can  find  the  dimension¬ 
less  components  of  stresses 


3*  ))  exp(/;£vl),  :  -  '.(C?— <-■)*  - 1  ;,(p))e':p,/ra'-.) 


where 


::*  =  v(C7— f  3,(p))  e\p  (//;  /!),  -,,=0,57  exp  ( 


//r-, 


(I. 15) 


-i(o)  =0..-. 


-  ml  *  j  dp 


ll.lG) 


( 


The  value  of  CJ  in  (1.16)  is  obtained  from  (1.11). 

Boundary  conditions  and  dimensionless  quantities  take  the  form 


37 


7'*(i)-nn=2.v(C,  rM(')).  /*(1)=0.  /*(:o)-«o 

(1.17) 

2.  At  a  small  value  of  parameter  v  we  set 


where 


r*=v7’.  Q*  =  -^"Q,  C, Q  (2.1) 


/=/o-r>/i  +  ''V2— 

T=7o+vr,  -rv*r3^  ...;  ^'  =/='o2-v7',tvV'’,+  ... 

go=  =^ou-f^<''oi+''%2  +  - 

g-2o(3/»^oV-i  /wf>'^7*)  =  Qo-t-vQi  f  v='Q,+  ... 

c, = J (^'  (>  -7)/ + (3+  = 

?U 

m—  1 

-1“ ‘>^C,2  "t"  ...*t  «>=^0 

Considering  relations  (2.2)  at  v  =  0,  (1.13)  and  (1.17)  give  the  system 
of  equations 


dfo  ^  :'?o --/u  .  ^-0  -Vo  -  pop. 

flfp  2'>  '  2t<* 

( 

^  1  — /n)p/-Vfo7'u-.3(5  -3/w)/,37',j-(3— /«)p-^^ru  -2p='/'^) 


dC'  Q„ 

and  boundary  conditions 


/■;(1)=0;  /■o(Po)=«o:  7'u(Po)«<o(P«)  =  2aiCio 


where 


7'j(lK|l)  =  »,(2C,.+  f 


+  ?■  ‘3»=2:-<3'i+ 


I  rii 

^^0  1%^  I  ‘“0  (3(  1  -  j)f -  (3  7  7)  f''Po)  ^P'  “’0  =  ^0 
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Solving  the  first  three  equations  of  (2.3)  and  satisfying  the  first  two 
boundary  conditions  of  (2.4),  we  find 


fo-C,(^-oy  F,=  2C, 


where  C3  =  PoUo/(l  -  pS)  and  is  an  arbitrary  constant. 

Substituting  (2.6)  in  the  last  equation  of  system  (2.3)  and  making 
necessary  simplifications,  we  get  a  nonlinear  differential  equation  for 

To 


(ci  (a  -  ,7)  +  W  ~  7 


whose  integral  takes  the  form 


ci(3+^)  +  T-(e;Tb)'- 


where  is  an  arbitrary  constant. 

Constants  and  Cg  are  determined  from  the  last  two  conditions  of  (2.4) 
which,  after  substitution  of  fo  and  (Po/  take  the  form 


1 


C  r  -  -  r  ^  3  ^  ^  do 

^  *  1-  J  V  o*)  T'o 


At  m  =  1,  we  find  from  (2.8) 


..  C--U)-) 


and  conditions  (2.9)  take  the  form 


C.  :  —  iC’s'.hl'V'o  +  ”  ^Q!’-2 


At  m  =  0,  we  get 


- T  _ 

To  ■■=  |/  ‘iry/ 472 _  ( Cj-^-Qo- )• 


and  at  m  =  0.5 


(•:.ii) 


(2.1:2) 


39 


I  (■-.'5TQ-7^4-K(cVrCV 


(34. I  ^.V) 


The  problem  has  been  numerically  studied  at  Pi  =  0.3,  Pa  =  0.2,  a  =  10, 
b  =  12.5,  Ui  =  20,  V  =  0.025  (Uo  =  1-6,  po  =  0.8)  for  different  values 
of  m:  m  =  0.1,  0.2,...,  0.9.  Values  of  P*/Co  have  been  calculated, 
where  P*  =  P/(kijb^),  P  is  the  force  of  pressing 


P  = 


j  H*V \  -  H*'-  '  -  "'o''"  (r't 

0 

=  -  =*(po>  'OA"'- 


A  “  '  2 

C'o-  ’) 


^u(r-'o)"“ot''o) 

(2.14) 


Fig.  2  shows  how  P*/i  depends  on  the  quantity  m. 

3.  In  the  case  of  external  introduction  (Fig.  3)  R*  =  1  -  vuq  exp  (vC) 
and  the  boundary  conditions  take  the  form 

U(b,  ^)=  -  Wij'exp  u(ii,  j)^0 

Formulas  (1.1)-(1.8),  (1. 11)“(1.16) ,  (2.1)-(2.3)  hold  in  this  case, 

while  the  remaining  fomulas  change 


T"(l)u< 

*(!)  = 

=  -Vr(6;- 

cd 

r(i)='-  --«o. 

=0 

ro(  i)=— i/o- 

^  IJ'  .  u 

)"'()(f'u) 

T„(  I  )o.o(  1 )  == 

3/o 

"o  - - 

' ! 

1 

(3.4) 

ro” 

/-o  =  2L', 

(3.5) 

/ 

V:., 

P  / 

a(-  ■ 

f-'/ 

Po  -  3r/' 

'(p{ 

3  ■  ' 

•  (1 

•i-  '  > 

.  L*  •  • 

)  r"'?) 

(3.6) 

(3.1) 

(3.2) 
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f  m  =  1,  then 

Cs+Qp*  (3.  in 

^  0 - 

0 

C5+C.pi]=-8p,C,,  -  2!t,C3(3+f'5)/Po  (3-12) 
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ELECTRICAL  TRANSIENT  PROCESSES  IN  THE  VAPOR- GAS  SHELL  UNDER  CONDITIONS  OF 
ELECTROLYTIC  ANODIC  HEATING  WITH  RECTANGULAR  VOLTAGE  PULSES 
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[Article  by  S.  I.  Galanin,  V.  I.  Ganchar,  and  E.  G.  Dmitriyev,  Institute  of 
Applied  Physics,  Moldavian  SSR  Academy  of  Sciences] 

[Text]  Anodic  electrolytic  heating  (AEH)  is  accompanied  by  film  boiling  of 
the  electrolyte  solution  on  the  surface  of  the  anode.  The  two-phase  system 
which  results  is  electrically  conductive,  and  the  existence  of  a  vapor  film 
(shell)  around  the  anode  results  from  the  joulian  heat  release  of  the  current 
flowing  through  it.  One  of  the  methods  of  obtaining  information  about  the 
nature  of  the  current  carriers  in  the  gaseous  phase  is  analysis  of  the 
dependence  of  the  electrical  conductivity  of  the  vapor-gas  shell  on  the 
applied  voltage.  However,  such  dependence  cannot  be  obtained  directly  from 
the  volt- ampere  characteristic  of  the  stationary  process  of  AEH,  for  in 
stationary  conditions  the  change  in  electrical  resistance  of  the  shell  under 
change  in  the  voltage  results  primarily  from  change  in  the  thickness  of  the 
shell  [1].  The  required  dependence  of  the  shell  conductivity  on  the  applied 
voltage  may  be  obtained  by  pulsed  measurement  methods.  In  the  most  simple 
experimental  layout  for  a  pulse  technique,  it  is  sufficient  to  measure  the 
time  constant  of  decay  of  the  discharge  current  of  a  capacitor  of  given 
capacitance  through  the  shell.  Such  a  simple  method  is  only  useful  when  the 
shell  represents  a  purely  active  resistance.  In  reality,  when  the  electric 
field  is  turned  on  in  pulses,  the  phase  boundaries  of  the  regions  near  the 
electrodes  are  charged  and  dual  electric  layers  are  created  on  them,  which  may 
be  approximated  as  capacitors  having  a  certain  overall  capacitance  Cj,.  In 
order  to  clarify  the  influence  of  this  capacitance  on  the  measurements  of 
electrical  conductance  of  the  shell,  it  is  necessary  to  know  its  magnitude. 
This  is  the  task  of  the  present  work. 

In  order  to  measure  the  magnitude  of  the  overall  surface  capacitance  Cn, 
rectangular  voltage  pulses  were  used.  A  steel  anode  was  heated  in  a 
stationary  AEH  process  up  to  T  ^  800® C,  and  then  the  voltage  was  switched  off. 
A  single  rectangular  voltage  pulse  was  applied  to  the  thus -prepared  current- 
free  shell  in  a  time  much  shorter  than  its  persistence  time.  The  results  of 
the  action  of  this  pulse  on  the  electrolytic  cell  with  and  without  an  anode 
vapor  shell  present  in  it  were  qualitatively  different  (see  Fig.  1):  the 
electrolytic  cells  does  not  distort  the  original  shape  of  the  imposed  pulse;  a 
distortion  is  only  observed  when  the  current- free  anode  shell  is  present  in 
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the  system.  This  experimental  fact  suggests  that  the  presence  of  the  current- 
free  anode  shell  in  the  electrolytic  cell  brings  about  an  additional 
capacitance  in  the  system,  sufficient  to  distort  significantly  the  shape  of 
the  imposed  pulse. 


Fig.  1.  Comparison  of  current  and  voltage 
oscillograms  obtained  with  a  rectangular 
voltage  pulse  acting  on  an  electrolytic 
cell  in  the  presence  (1)  and  absence  (2) 
of  a  current- free  vapor  shell  in  it. 

Key:  a.  /is. 


It  is  easy  to  support  this  conclusion  with  calculation.  For  this,  we  repre¬ 
sent  the  electrolytic  cell  with  current- free  anodic  shell  present  in  it  as  an 
equivalent  substitution  circuit  (Fig.  2).  The  current  and  voltage  measured  by 
the  oscillograph  may  be  obtained  from  the  system  of  equations: 


C  0 

=  /«(0)  =  0, 


(1) 


where  Uq  is  the  potential  of  the  source  of  rectangular  pulses;  lo  is  the 
current  flowing  through  the  electrolytic  cell,  Ir  and  Ip  are  the  currents 
flowing  through  the  active  resistance  of  the  phase  interface  surfaces  Rn  and 
the  overall  surface  capacitance  respectively;  R^,  Rj,  and  are  the 

internal  resistance  of  the  source  and  the  active  resistances  of  the  electro¬ 
lyte  and  shell,  respectively.  From  (1),  we  obtain  for  the  sought  quantities 
U(t)  and  Io(t): 
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where 


=  /?2  =  ^3  +  ^o6  +  ^h,  Rz  —  Ri+RoO  +  Ru-^ 

-f/?n. 


Calculation  by  (2)  and  (3)  yields  time  functions  U(t)  and  Io(t),  similar  to 
the  experimental.  The  observed  distortion  in  shape  of  the  imposed  pulse  is 
achieved  when  Cn  ^  2000  fiF  for  values  of  Rj,  Rj,  and  R3  on  the  order  of 
0.1-1  ohms  and  duration  of  the  imposed  pulse  r  =  500  /is. 


Fig.  2.  Equivalent  electrical  substitution 
circuit  of  an  electrolytic  cell  with 
current-free  vapor  shell  present  in  it. 


The  surface  capacitance  of  the  phase  interface  boundary  can  also  be  evaluated 
from  other  considerations.  As  an  example,  consider  an  electrolyte-vapor 
interface  surface.  Under  the  action  of  the  imposed  voltage,  the  phase 
boundary  is  charged  until  the  electric  field  of  the  surplus  charge  accximulat- 
ing  on  the  surface  cancels  out  the  external  field.  The  characteristic 
dimension  of  the  charge  separation  will  be  the  Debye  screening  radius,  which 
is  connected  to  the  equilibrium  charge  density  within  the  electrolyte  n  by  the 
familiar  relation: 


'‘x)==(e7'/4  ;te%)V2^ 

Here,  T  is  the  temperature  of  the  electrolyte  near  the  surface,  expressed  in 
energy  units;  e  is  the  charge  of  the  electron;  e  is  the  dielectric 
permittivity  of  the  electrolyte  near  the  phase  boundary.  The  surface  layer  of 
electrolyte,  r^  thick,  may  be  approximated  as  a  flat  capacitor,  the  capaci¬ 
tance  of  which  is : 


Cfi —  Dt  (  5  ) 

where  Cq  ^  8. 85- 10'*^  F/m  is  the  dielectric  constant.  As  a  result,  we  find  that 
Cj,  ~  yn,  and  as  the  electrolyte  concentration  varies  from  13  to  38  percent  the 
specific  surface  conductivity  of  the  electrolyte -vapor  phase  boundary  changes 
from  0.82  to  1.5  F/m^.  In  the  conditions  under  which  the  measurements  were 
taken,  the  anode  surface  was  2.7-10*^  m^;  thus,  the  total  capacitance  of  the 
phase  surface  should  change  from  2200  to  3900  fiF  as  the  electrolyte 
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concentration  varies  in  the  indicated  range.  Such  dependence  of  the  surface 
capacitance  on  the  electrolyte  concentration  is  experimentally  confirmed:  as 
the  electrolyte  concentration  rises,  the  shape  of  the  imposed  voltage  pulses 
is  more  heavily  distorted  (Fig.  3). 


Fig.  3.  Oscillograms  of  voltage  at  various  concentra¬ 
tions  of  electrolyte  NH4NO3:  1-32  percent,  2  - 
13  percent.  Solid  lines  -  calculation,  points  - 
measurement  results . 

Key: 

a.  V  b.  , 


As  for  the  nature  of  the  influence  of  the  surface  capacitance  of  the  phase 
boundaries  on  the  pulse  measurement  results,  it  may  be  variable.  For  example, 
when  low-voltage  pulses  are  applied  to  the  current-free  anode  shell  and  the 
potential  Uq  is  less  than  a  certain  critical  value  U„,  there  are  no  conduction 
currents  through  the  vapor  shell  [2],  and  the  surface  capacitance  Cn  in  the 
equivalent  substitution  circuit  is  connected  in  series  with  the  proper 
capacitance  of  the  shell  Since  Cq(J«  0^,  the  overall  capacitance  of  the 

system  is  C  ==  Thus,  the  influence  of  the  surface  capacitance  is 

eliminated,  making  it  possible  to  measure  the  proper  capacitance  of  the  anode 
shell.  Such  measurements  were  performed  in  [3];  according  to  the  data  of  this 
work,  (150-300)  pF.  A  different  situation  occurs  when  measuring  the 

conductance  of  the  anode  shell:  when  a  capacitor  charged  to  potential  Uq  > 
is  connected  to  the  shell,  a  discharge  current  flows  through  it;  the 
capacitance  of  the  phase  interface  surface  is  added  to  the  capacitance  being 
discharged.  The  capacitance  of  the  phase  interface  surface  is  on  the  order  of 
thousands  of  microfarads,  while  the  electrolyte -vapor  phase  surface  is  able  to 
discharge  in  a  time  on  the  order  of  hundreds  of  microseconds.  Consequently, 
in  order  to  eliminate  the  influence  of  the  capacitance  of  the  phase  boundary 
on  the  measurements,  the  discharging  capacitance  must  be  chosen  in  the  range 
of  values  much  larger  than  the  capacitance  of  the  phase  boundary,  and  the 
current  decay  curve  must  be  registered  in  the  millisecond  time  range  in  a 
period  longer  than  the  charge  time  of  the  phase  interface  surface. 
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MECHANICS  OF  GASES,  LIQUIDS,  SOLIDS 


INVESTIGATION  OF  AN  UNDERWATER  ELECTRIC  EXPLOSION  INITIATED  BY  FERROMAGNETIC 
PARTICLES.  I.  PHYSICAL  MODELING 

907F0366B  Kishinev  ELEKTRONNAYA  OBRABOTKA  MATERIALOV  in  Russian  No  2,  1990 
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[Article  by  L.  Z.  Boguslavskiy ,  L.  M.  Bondarets,  V.  V.  Kucherenko,  and  V.  V. 
Shamko,  Electrohydraulic  Design  Development  Office,  Ukrainian  SSR  Academy  of 
Sciences] 

[Text]  The  results  of  experimental  studies  of  an  electrical  explosion  in 
ferromagnetic  suspensions  [1,  2]  show. that,  depending  on  the  conduction  and 
concentration  of  the  solid  phase,  both  the  conditions  of  formation  of  the 
channel  and  the  discharge  regime  vary.  This  should  result  in  change  in  the 
gas  dynamic  characteristics  of  the  underwater  electrical  explosion  (UEE)  and, 
accordingly,  the  effectiveness  of  transformation  of  electrical  energy  into 
mechanical  [3,  4].  Therefore,  the  present  work  is  devoted  to  developing  the 
methods  of  calculation  of  the  energy  and  gas  dynamic  characteristics  of  an 
electric  explosion  in  a  liquid  containing  ferromagnetic  inclusions. 

In  the  first  part  of  the  work,  the  methods  of  physical  modeling  are  used  to 
investigate  the  electrical,  energy,  and  kinematic  characteristics  of  the  UEE 
with  initiation  by  ferromagnetic  particles,  localized  in  lengthwise  bridges. 

The  experiments  were  carried  out  in  a  tank  (0.9  x  0.9  x  0.6  m*)  with  portholes 
of  organic  glass,  filled  with  distilled  water  of  electric  conductance 
CTo  =  2.5«10'^  O'  X  m’’  (the  layout  is  similar  to  that  of  [5]).  In  order  to 
reduce  the  cumulative  effect  of  the  electric  explosion  cavity  (EEC) ,  the 
negative  electrode  had  a  rectangular  opening  of  0.03  x  0.23  x  0.1  m^  (the 
external  appearance  of  the  electrode  system  with  ferromagnetic  bridge  is  shown 
in  Fig.  1).  The  distance  between  electrodes  in  all  cases  was  1  =  40  mm,  and 
the  discharge  occurred  at  depth  of  0.3  m. 

The  bridges  were  created  from  ferromagnetic  powder  SF-35  with  initial  electric 
conductivity  Oq  =  10'^  fl'»m''  as  follows.  In  a  separate  vessel,  a  10  percent 
suspension  of  powder  SF-35  was  prepared  from  distilled  water,  and  a  capillary 
was  used  to  supply  this  to  the  interelectrode  gap  (lEG)  with  a  lengthwise 
magnetic  field,  created  by  a  special  coil  [1] .  When  the  suspension  enters  the 
region  of  the  field,  the  solid  particles  line  up  along  the  field  lines  and 
form  a  bridge  spanning  the  entire  lEG  (see  Fig.  1). 
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Fig.  1.  External  appearance  of  electrode 
system  and  ferromagnetic  bridge:  1  - 
negative  electrode;  2  -  positive 
electrode;  3  -  ferromagnetic  bridge. 


The  electrotechnical  parameters  (charge  voltage  of  the  capacitor  bank  Uq,  its 
capacitance  C,  the  inductance  of  the  discharge  circuit  L)  are  presented  in 
Table  1.  For  each  particular  set  of  parameters  Uq,  C,  L,  at  least 
13-15  discharges  were  performed,  and  the  obtained  results  were  averaged. 


Table  1.  Electrotechnical  parameters  of  the  discharge  circuit 
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Key:  a.  Number  of  discharge 

b.  V  c.  F  d.  J 

NOTE:  L,  lO"^  =  3.5  H. 

The  studies  were  done  by  the  method  simultaneous  registration  of  electric  and 
space- time  characteristics  of  the  discharge  channel  and  the  EEC.  The  current 
was  registered  by  a  low-resistance  coaxial  shunt,  the  voltage  drop  across  the 
discharge  gap  by  a  capacitive  divider,  the  signals  from  which  were  sent  to  the 
S8-17  oscillograph.  The  inductive  component  of  the  voltage  drop  was  compen¬ 
sated.  The  discharge  channel  was  photographed  by  the  high  speed  camera  VFU-1, 
acting  as  a  photoregister  (mirror  rotational  speed  30,000  rpm) ,  the  cavity  was 
photographed  by  the  SFR-2M  with  an  adapter  reducing  the  mirror  turning  speed 
by  two  orders  of  magnitude.  Highlighting  during  investigation  of  the  channel 
stage  was  provided  by  the  pulse  lamp  IFK-2000,  and  after  the  discharge  the 
constant  light  source  KG- 110  was  used. 


The  current  and  voltage  oscillograms  (similar  to  those  presented  in  [1]  were 
used  to  determine:  in  the  channel  formation  stage,  the  breakdown  delay  time 
tj,  the  decrease  in  magnitude  of  the  initial  potential  AUq,  the  energy  losses 
AW  (averaged  values  of  these  are  given  in  Table  2  for  the  different  discharge 
regimes),  and  in  the  active  stage  the  energy  generated  in  the  channel  W(t), 
the  time  for  its  generation  r,  and  the  electric  power  N(t) .  Deviations  of  the 
investigated  quantities  from  the  mean  values  usually  did  not  exceed  30  percent 
for  tj,  50  percent  for  AU,,  10  percent  for  r,  and  15  percent  for  AW/Wo. 


Table  2.  Averaged  time  and  energy  characteristics 
of  discharges  1-3 
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A  characteristic  of  the  formation  of  the  channel  when  the  discharge  is 
initiated  by  powder  SF-35  is  the  fact  that,  within  20-30  fis  after  voltage  is 
applied  to  the  lEG,  several  luminous  regions  of  size  1-2  mm  are  created  at  its 
center,  the  number  of  which  increases  in  the  lengthwise  direction  over  the 
time  (Fig.  2a,  b) .  Occasionally,  two  luminous  chains  appear, 
spanning  the  entire  gap  at  a  certain  instant. 

The  mechanism  of  formation  of  the  discharge  channel  (breakdown  stage)  when  the 
discharge  is  initiated  by  conductive  ferromagnetic  inclusions  has  been 
discussed  in  [1,  2,  6]  and  is  explained  by  the  redistribution  of  the  electric 
field  in  the  ga.p  as  a  result  of  the  presence  of  the  bridge,  resulting  in 
concentration  of  current  and  joulian  heat  release  in  individual  regions  on  the 
surface  of  the  inhomogeneity.  The  latter,  in  the  opinion  of  [1,  2,  6],  may 
result  in  initiation  of  a  discharge  in  these  zones.  It  may  be  that  a  similar 
breakdown  mechanism  occurs  in  a  suspension  based  on  the  powder  SF-35,  since 
its  electric  conduction  is  greater  than  that  of  water.  However,  other 
mechanisms  may  be  at  play  (breakdown  along  the  surface  of  the  bridge,  break¬ 
down  of  gas  bubbles  ever  present  in  the  liquid,  etc.). 

It  is  difficult  to  distinguish  in  time  and  space  the  stages  of  formation  of 
the  discharge  channel  (prior  to  breakdown)  and  the  channel  proper,  as  is 
obvious  from  Fig.  2a.  Toward  this  end,  the  discharge  channel  was  photographed 
without  lighting  (see  Fig.  2b)  and  the  appearance  of  an  intense  diffuse  light 
around  the  channel  in  the  second  frame  was  coordinated  (on  the  analogy  of  a 
discharge  without  initiation)  with  the  start  of  a  sharp  rise  in  current  on  the 
oscillograms.  Thus,  the  electric  and  kinematic  characteristics  of  the 
discharge  when  initiated  by  ferromagnetic  bridge  were  synchronized. 
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Fig.  2.  Motion  picture  records  of  the  channel  and  cavity  for  the 
first  discharge:  a,  b  -  discharge  channel  with  UEE  initiated  by 
powder  SF-35;  c  -  cavity  with  UEE  initiated  by  copper  wire, 
diameter  0.1  iran. 

Key:  d.  /is  e.  ms 


The  resulting  discharge  channel  is  not  homogeneous  either  along  the  lengthwise 
or  radial  coordinate  (its  glow  is  not  uniform),  making  it  hard  to  determine 
its  kinematic  characteristics .  Figure  3  shows  the  functions  of  the  radius  of 
the  channel  a(t)  and  the  speed  of  its  expansion  a(t)  for  the  first  discharge. 
The  time  is  counted  off  from  the  instant  of  appearance  of  glow  in  the  lEG  (the 
stage  prior  to  breakdown  is  also  partly  included) .  These  curves  are  basically 
qualitative  in  nature,  due  to  the  difficulties  of  exact  determination  of  the 
radius  of  the  channel  (in  the  present  case,  it  is  determined  as  the  radius  of 
the  brightest  portion  of  the  channel  for  a  particular  cross  section  that  is 
sufficiently  visible  in  the  motion  picture  records,  see  Fig.  2a).  As  is 
evident  from  Fig.  3,  the  parameters  of  the  channel  at  the  instant  of  closing 
of  the  plasma  branch  of  the  lEG  (t  ~  30  /is),  i.e.,  the  initial  conditions  for 
the  channel  stage,  amount  to  ^  1  mm,  ^  «  10  m/s. 

Comparison  of  the  experimental  curves  of  the  power  N(t)  ,  obtained  in  the 
present  work  for  electric  explosion  initiated  by  SF-35  powder  (Fig.  4, 
curve  2)  and  initiated  by  microconductor  or  high-voltage  liquid  breakdown  [3], 
revealed  a  qualitative  analogy.  Previously  [3],  an  analogy  was  found  for  the 
electric  and  kinematic  characteristics  of  UEE  in  terms  of  the  parameter: 


characterizing  the  discharge  regime  and  representing  the  portion  of  energy 
released  during  the  first  half -period  of  the  current.  Here,  A  is  the  spark 
constant,  which  for  an  electric  explosion  initiated  by  microconductor  is 
A  =  0.25#10^  V^ts/m^. 


Fig.  3.  Time  curves  of  radius  of  channel 
a(t)  and  speed  of  its  expansion  3.(t)  for 
the  first  discharge:  1  -  experiment,  2  - 
calculation. 

Key:  a.  m/s  b.  /is 


In  the  present  work,  the  parameter  ?y  was  determined  experimentally  as  the 
ratio  between  the  energy  Wj  released  during  the  first  half -period  of  the 
current  and  the  stored  energy  Wq  =  It  turned  out  that,  for  unchanged 

initial  electrotechnical  parameters  of  the  channel  and  physico -mechanical 
properties  of  the  working  medium,  the  parameter  changes  from  one  discharge  to 
another  (for  an  UEE  initiated  by  microconductor,  the  regime  is  uniquely 
characterized).  The  range  of  change  in  ly  for  discharges  1-3  is  shown  in 
Table  2.  While  for  conducting  powders  [1,  2,  6]  this  was  due  to  change  in 
concentration  of  the  solid  phase  in  the  ferromagnetic  suspension  and,  thus, 
the  initial  diameter  of  the  bridge,  for  the  nonconducting  powder  SF-35  this 
may  be  explained  by  an  increase  in  length  of  the  channel  as  compared  to  the 
length  of  the  lEG  (see  the  above  expression  for  »/).  Values  of  the  equivalent 
channel  lengths  Ij,  corresponding  to  the  UEE  initiated  by  microconductor  and 
calculated  for  different  discharge  regimes,  are  presented  in  Table  2.  The 
above-given  results  of  experimental  investigations  of  the  energy  character¬ 
istics  of  an  electric  explosion  initiated  by  SF-35  powder  show  that  they  can 
be  determined  in  accordance  with  the  technique  [4]  developed  for  an  ordinary 
underwater  spark  discharge.  However,  in  calculating  the  parameter  i;,  instead 
of  the  length  of  the  gap  1  it  is  necessary  to  use  21  <  Ij  <  31. 
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Fig.  4.  Power  curves  for  the  first 
discharge:  1  -  experiment,  2  -  calcula¬ 

tion  by  [4] . 

Key:  a.  W  b.  /iS 


We  investigated  the  influence  of  the  opening  in  the  negative  electrode  on  the 
kinematics  of  the  cavity.  Since  a  poor  quality  of  photograph  of  the  EEC  is 
obtained  when  the  UEE  is  initiated  by  powder  SF-35,  Fig.  2c  presents  photo¬ 
graphs  for  the  first  discharge  initiated  by  copper  conductor  of  diameter 
0.1  mm.  During  the  first  pulsation,  the  cavity  is  dome  shaped,  and  becomes 
slightly  deformed  in  the  closure  stage.  It  has  been  established  that  the 
shape  of  the  cavity,  its  maximum  volume,  and  the  period  of  the  first  pulsation 
are  almost  unaffected  by  the  opening,  which  basically  influences  the  speed  of 
movement  of  the  EEC  boundaries  (the  cavity  collapses  more  intensely  if  the 
opening  is  missing,  see  Fig.  5a,  b)  and  the  number  of  compression  pulses 
generated  by  the  cavity.  Thus,  when  an  opening  is  present,  the  EEC  produces  a 
single  compression  wave,  and  when  it  is  absent  (closed  by  a  metal  plate),  it 
produces  two  waves.  Analysis  of  the  kinematic  characteristics  of  the  EEC 
(height  of  the  dome  djL(t)  ,  its  speed  djjit)  in  the  perpendicular  direction  (a) 
and  the  radius  of  the  cavity  a^(t),  the  speed  aj|(t)  in  the  direction  parallel 
to  the  negative  electrode.  Fig.  5)  reveals  that  the  functions  more  sensitive 
to  presence  of  an  opening  are  dj^(t),  dj^(t) .  The  kinematic  characteristics  of 
the  EEC  formed  during  UEE  initiated  by  powder  SF-35  for  the  third  discharge 
are  shown  in  Fig.  5c. 

In  the  concluding  stages  of  the  closure  (see  Fig.  2c),  the  cavity  moves  in  the 
direction  of  the  free  surface  and  collapses  to  form  a  cuimulation  stream  of 
liquid,  4  cm  above  the  point  of  contact  between  the  channel  and  the  plate 
covering  the  opening.  At  the  point  of  impact  of  the  stream  on  the  plate,  a 
conical  funnel  is  formed  with  radius  of  depression  at  the  contact  point 
1-2  mm.  Because  the  regions  of  influence  of  the  channel  and  the  cavity 
(stream)  are  at  a  distance  from  each  other,  it  was  possible  to  separate  their 
mechanical  action  on  the  plate.  It  was  found  that  the  bending  of  the  plate 
(using  copper  and  steel  strips,  6  =  0.25-1  mm  thick)  at  the  point  of  action  of 
the  shaped  stream  of  liquid  was  several  millimeters  larger  than  at  the  point 
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of  contact  of  the  channel.  Thus,  for  a  copper  plate  with  5  -  0.25  mm,  the 
bending  was  respectively  10  and  7  mm.  This  result  is  consistent  with  the  data 
of  [7,  8]. 


Fig.  5.  Kinematic  characteristics  of  cavities:  a,  b  -  initiated 
by  copper  conductor,  0.1  mm  diameter,  first  discharge;  c  ~ 
initiated  by  powder  SF-35,  third  discharge;  1,2-  experiment 

(primed  ntunbers  -  radius),  3  -  calculation; - electrode  without 

opening . 


Thus,  the  initiation  of  a  discharge  by  nonconducting  ferromagnetic  particles 
satisfactorily  stabilizes  its  energy,  kinematic,  and  time  characteristics. 


Bibliography 

1.  Demina,  V.  M.  ,  Lyubimov,  A.  D.  ,  Satanovskiy,  V.  V.,  "Experimental 
Studies  of  the  Breakdown  of  a  Ferromagnetic  Suspension  in  an  External 
Magnetic  Field,"  ELEKTRONNAYA  OBRABOTKA  MATERIALOV,  No  4,  1985, 

pp.  52-56. 

2.  Nazarova,  L.  L. ,  Malikov,  P.  D.,  Reztsov,  V.  F. ,  "On  Transient  Processes 
in  Electric  Discharge  in  Liquid  Containing  Ferromagnetic  Inclusions  in 
an  External  Lengthwise  Magnetic  Field,”  "Fizicheskiye  osnovy 
elektricheskogo  vzrjrva:  Sb.  nauch.  tr."  [Physical  Foundations  of  an 
Electric  Explosion:  Collection  of  Scientific  Papers],  Kiev,  1983, 

pp.  60-64. 

3.  Krivitskiy,  Ye.  V.,  Shamko,  V.  V.,  "Perekhodnyye  protsessy  pri 
vysokovoltnom  razryade  v  vode”  [Transient  Processes  in  High-Voltage 
Discharge  in  Water] ,  Kiev,  1979. 

4.  Kucherenko,  V.  V. ,  "Nestatsionarnyye  gazodinamicheskiye  techeniya  pri 
podvodnom  elektrovzryve"  [Nonstationary  Gas  Dynamic  Currents  in  an 


54 


Underwater  Electric  Explosion],  Author's  Abstract  of  Dissertation  for 
Candidate  of  Technical  Sciences,  Kiev,  1988. 

5.  Shamko,  V.  V.,  Vovchenko,  A.  I.,  "Influence  of  Boundary  Surfaces  on 

Development  of  the  Vapor-Gas  Cavity  in  an  Underwater  Spark  Discharge," 
"Gidromekhanika:  Sb.  nauchn.  tr.”  [Hydromechanics:  Collection  of 

Scientific  Papers],  Issue  34,  Kiev,  1976,  pp.  52-60. 

6.  Malikov,  P.  D.  ,  Nazarova,  L.  L.  ,  Reztsov,  V.  F.  ,  "A  Qualitative  Model  of 
Initiation  of  an  Electric  Discharge  in  a  Conductive  Liquid  by 
Inhomogeneities,"  ZHTF,  Vol  55,  No  3,  1985,  pp.  616-618. 

7.  Khasiguti  Kiyekhito,  Yamatoto  Takekhisa,  "Movement  of  a  Gas  Bubble  and 
Impact  Force  During  Discharge  of  a  Capacitor  in  Liquid,"  KIKAY  TO  KOGU, 
Vol  13,  No  11,  1969,  pp.  146-154. 

8.  Chachin,  V.  N. ,  "Elektrogidravlicheskaya  obrabotka  mashinostroyitelnykh 
materialov"  [Electrohydraulic  Machining  of  Materials],  Minsk,  1978. 

COPYRIGHT:  Izdatelstvo  "Shtiintsa,"  "Elektronnaya  obrabotka  materialov,"  1990 

-  END  - 


55 


This  is  a  U.S.  Government  publication.  Its  contents  in  no  way  represent  the 
policies,  views,  or  attitudes  of  the  U.S.  Government.  Users  of  this  publication  may 
cite  FBIS  or  JPRS  provided  they  do  so  in  a  manner  clearly  identifying  them  as  the 
secondary  source. 


Foreign  Broadcast  Information  Service  (FBIS)  and  Joint  Publications  Research  Service  (JPRS) 
publications  contain  political,  military,  economic,  environmental,  and  sociological  news,  commentary, 
and  other  information,  as  well  as  scientific  and  technical  data  and  reports.  All  information  has  been 
obtained  from  foreign  radio  and  television  broadcasts,  news  agency  transmissions,  newspapers,  books, 
and  periodicals.  Items  generally  are  processed  from  the  first  or  best  available  sources.  It  should  not  be 
inferred  that  they  have  been  disseminated  only  in  the  medium,  in  the  language,  or  to  the  area  indicated. 
Items  from  foreign  language  sources  are  translated;  those  from  English-language  sources  are  tran¬ 
scribed.  Except  for  excluding  certain  diacritics,  FBIS  renders  personal  names  and  place-names  in  accor¬ 
dance  with  the  romanization  systems  approved  for  U.S.  Government  publications  by  the  U.S.  Board 
of  Geographic  Names. 

Headlines,  editorial  reports,  and  material  enclosed  in  brackets  [  ]  are  supplied  by  FBIS/JPRS. 
Processing  indicators  such  as  [Text]  or  [Excerpts]  in  the  first  line  of  each  item  indicate  how  the 
information  was  processed  from  the  original.  Unfamiliar  names  rendered  phonetically  are  enclosed  in 
parentheses.  Words  or  names  preceded  by  a  question  mark  and  enclosed  in  parentheses  were  not  clear 
from  the  original  source  but  have  been  supplied  as  appropriate  to  the  context.  Other  unattributed 
parenthetical  notes  within  the  body  of  an  item  originate  with  the  source.  Times  within  items  are  as  given 
by  the  source.  Passages  in  boldface  or  italics  are  as  published. 


The  FBIS  DAILY  REPORT  contains  current  news 
and  information  and  is  published  Monday  through 
Friday  in  eight  volumes:  China,  East  Europe,  Soviet 
Union,  East  Asia,  Near  East  &  South  Asia,  Sub- 
Saharan' Africa,  Latin  America,  and  West  Europe. 
Supplements  to  the  DAILY  REPORTS  may  also  be 
available  periodically  and  will  be  distributed  to  regular 
DAILY  REPORT  subscribers.  JPRS  publications,  which 
include  approximately  50  regional,  worldwide,  and 
topical  reports,  generally  contain  less  time-sensitive 
information  and  are  published  periodically. 

Current  DAILY  REPORTS  and  JPRS  publications  are 
listed  in  Government  Reports  Announcements  Issued 
semimonthly  by  the  National  Technical  Information 
Service  (NTIS),  5285  Port  Royal  Road,  Springfield, 
Virginia  22161  and  the  Monthly  Catalog  of  U.S.  Gov¬ 
ernment  Publications  issued  by  the  Superintendent  of 
Documents,  U.S.  Government  Printing  Office,  Wash¬ 
ington,  D.C.  20402. 

The  public  may  subscribe  to  either  hardcover  or 
microfiche  versions  of  the  DAILY  REPORTS  and  JPRS 
publications  through  NTIS  at  the  above  address  or  by 
calling  (703)  487-4630.  Subscription  rates  will  be 


expect  a  30-day  delay  in  receipt  of  the  first  issue. 

U.S.  Government  offices  may  obtain  subscrip¬ 
tions  to  the  DAILY  REPORTS  or  JPRS  publications 
(hardcover  or  microfiche)  at  no  charge  through  their 
sponsoring  organizations.  For  additional  information 
or  assistance,  call  FBIS,  (202)  338-6735, or  write 
to  P.O.  Box  2604,  Washington,  D.C.  20013. 
Department  of  Defense  consumers  are  required  to 
submit  requests  through  appropriate  command  val¬ 
idation  channels  to  DIA,  RTS-2C,  Washington,  D.C. 
20301.  (Telephone;  (202)  373-3771,  Autovon; 
243-3771.) 

Back  Issues  or  single  copies  of  the  DAILY 
REPORTS  and  JPRS  publications  are  not  available. 
Both  the  DAILY  REPORTS  and  the  JPRS  publications 
are  on  file  for  public  reference  at  the  Library  of 
Congress  and  at  many  Federal  Depository  Libraries. 
Reference  copies  may  also  be  seen  at  many  public 
and  university  libraries  throughout  the  United 
States. 


SUBSCRIPTION/PROCUREMENT  INFORMATION 

provided  by  NTIS  upon  request.  Subscriptions  are 
available  outside  the  United  States  from  NTIS  or 
appointed  foreign  dealers.  New  subscribers  should 


